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1.0 INTRODUCTION AND BACKGROUND 


Hydrologists and water resource planners are continually 
faced with developing two types of information concerning 
watersheds: 1) Management information which deals with the 

response of the watershed to specific individual event (s), 
and 2) Planning information necessary to define the expected- 
response of the watershed to statistically recurring peak 
events. Management information is required for the real-time 
optimization of supply versus demand; planning information is 
needed for the optimal sizing of waterworks. 

Typical- -planning events treated by water resources agencies 
Include: probability of recurrence of peak flow rates (flood 

flows), in order to size spillways and determine .floodways: 
recurrence -duration of the high flow event (volume of water), 
in order to determine storage required to prevent flooding; 
recurrence-duration of the low flow event to maintain depend- 
able basin yield; and other water-related characteristics such 
as sediment yield and water quality. 

Due to the pressures of urbanization, one of the most important 
of the planning events above as regards both manpower employed 
and value resulting from optimal operations,, is the prediction 
of the peak flow event- flood frequency. 
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Nou ali^of the v/atershed phS^slcal- characteristics are 
equally important in determining the recurrence and magnitude 
of peak flow events . ‘It is thus fruitful to investigate in 
detail the following elements : 

1. Which of the watershed physical parameters are the 
most significant contributors to flood events; 

2. What are the errors in prediction committed by 
neglecting the less important parameters; 

/ 

3. What are the errors induced by imperfect knowledge 
of the principal parameters (sensitivity); 

and 

i|. How do the important "driver” parameters vary as a 
function of regional conditions. Understanding and 
- reduction to engineering practice of this last item 
would yield the important result of being able to 
specify the functional form, and hopefully the co- 
efficients, of a planning model for each region, 
following a procedure valid for all regions 
(regionalization) . 

The principal, conventional flood frequency planning models 
employ either 1) statistical extrapolation, particularly' for 
small-medium watersheds, or 2) digital parametric simulation 
modeling for large watersheds having high economic value. 

Both Of these methods require extensive historical rainfall 
and runoff records. Many watersheds do not have the required 
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length of records, particularly if the watershed is under- 
going rapid changes du“e' to urbanization; in this case, in 
fact, length of record is of very limited value. 


Another important question pertaining to planning models 
is being asked more .and more frequently by local .users.: 
specifically, what is the effect upon the hydrologic regime 
of planned modifications to ..the watershed, such as construction, 
deforestation, and reforestation, and how will these changes 
affect the size, type and location of .waterworks? . 

Rapid, repetitive survey techniques are required which relate 
the- physical characteristics of the watershed, particularly 
the surface characteristics, to information necessary for 
optimal planning of watershed development . If this informa- 
tion can be practically incorporated into planning models 
sensitive to the spatially distributed characteristics of the 
watershed (e,g,,.-y.egetatiye cp.ver,. impermeable/permeable areas, 
surface water, drainage pattern, and evapotranspiration po- 
tentlal), it will prove of major interest to: 

® Federal Agencies operating in the Water Resources 

field, such as the Corps of Engineers, the U.S, 

* * •; * * . ' > 

i . ***** 

Geological Survey, the USDA-Soil Conservation 

Service, the Agricultural Research Service, and the 
National Oceanographic and Atmospheric Administration. 
Federal Agencies whose mission is to guide and foster 
Water Resources Research, such as DOI - Office of 


O 



•Water Research anh Technolo.^y and the National JScience 
Fdundation' Intfernationai Hydrologic Program. 

® The developing nations wherein water resources develop- 
ment is a major influence in well-being and where in 
general the paucity of records, and the lack of a hy- 
drologic infra-structure* place critical importance 
on this effort. 

® Foreign governmental and extra-governmental agencies 
with water resources development orientation, such as 
the United Nations, the PAO, and the International Water 
Resources Association. 

1.1 State of the Art of Planning Models 

The historical evolution of planning models has progressed 
along three lines, leading to the formulation of three broad 
categories of methods, 

1.1.1 Method 1: The Empirical Approach 

The earliest to be devised, and perhaps still the most widely 
en^loyed, particularly in ungaged areas, are empirical formu- 
lations of the general type: 
q s A-^f(A)h(W)g(T) 

Where: 

q * peak flow rate from unit watershed area 
A = total watershed area 
f(A)* functional adjustment for area 



- 5 - 


h(W) « functional relationship of watershed physical 
characteristics, e.g. , slope, vegetative cover 

.g(T) = recurrence period, years 

n * empirically derived coefficient 

There exist perhaps a hundred models of this type in use in 
the U.S. and the world, particularly in ungaged areas. The 
significant models in general use in the world are presented 
in Table 1.* Their performances are compared graphically in 
Figures 1 and 2. As can be seen, there is a variation in the 
predicted peak discharge between models of at least one order 
of magnitude for the small and medium watersheds, increasing 
to two orders of magnitude for the large watersheds . 

This variation, per se, would not constitute a major problem 
if a unique model could be attributed to each region. In 
reality, even within the same geographic area, the model 
formulations vary significantly, as shown by Figure 3. An 
exan^le of the variability of the results of typical empirical 
regression models is given in Figure ^ for a large watershed 
in Maryland. There is a significant uncertainty in predictions 
for long recurrence periods which are usually related to de- 
sign of civil structures. Improvement in model accuracies in 
this region - and in many similarly developed regions - are 
of significant economic value. 

*lSource: Gray, Introduction to Hydrology 
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1.1,2 Method 2: The Statistical Approach 

If there is a sufficiently long historical record available.-, 

the random character of the flow peak can be treated, as a 

stochastic variable. There is an underlying assumption that 

the controlling random process is stationary, e.g., the 

character of the event remains unchanged by long term trends 

or other effects. Using daily records, various statistical 

parameters such as the mean, the variance, the coefficient of 

variability, and the skewness can be computed. By assuming a 

1 

given probability distribution such as Log Pearson III, Gurabel, 
or log-normal, the probability of the recurrence of a peak e- 
vent equal to or greater in magnitude than an arbitrary value 
can be determined. 

There are several fundamental problems involved with the pure- 
ly statistical approach. First, there is the obvious require- 
ment for extensive records in order to achieve a given level 
of confidence. This dependence is shown in Figure 5, which 
indicates the sensitivity of the record length in predicting 
the fifty year event for a particular region in the U.S. 
Percent of error due to this cause alone decreases from U0% 
to 1635 as the length of record increases from- 2 to 20 years 
a.t'90% confidence level. This assumes that the watershed has 
not changed significantly during the period of record. 
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AdditionaXly, there is considerable variation among pre- 
dictions from the principal statistical models in use and 
between the model results and actual records . Figure 6 is 
a comparison between the two principal distributions used 
in the U.S. and shows the magnitude of the error for two 
typical watersheds. 

In summary, the statistical method is not applicable to un- 
gaged watersheds, and is subject to potentially significant 
errors in watersheds undergoing change. 

1.1.3 Method 3: Semi-Empirical Macro Models - Rational Approach 

A useful and logical extension to these wholly empirical models 
is the addition of the rainfall component Iq, This component 

is the average rate of rainfall ”i" (nv/hr) which is observed 
to occur for a duration "d** (hrs) at a recurrence interval of 
"T" (years). Since a large number of rainfall records are 
generally available, it is reasonable to assume that the in- 
clusion of rainfall data should reduce the variance of the 
results. Semi -empirical methods are useful, although far 
from precise, where the coefficients relating runoff to rain- 
fall are well-known; however, extrapolation, of the coeffici- 
ents on a regional basis is 'difficult and the results are 
subject to considerable uncertainty. 

Rational or semi -empirical models have the general form: ‘ 
q = A«iro(T)g(T)h(W) 




Where: 


• q » flow per unit watershed area 
A = total watershed area’ 

i s rainfall rate of /a given -duration occurring 
at a given recurrence period T 

W = watershed parameters, e.g. Vegetative cover, 

• drainage density, slope, etc. 

g,h= functional relationships 

- * 1 - • 

n,m= empirically derived coefficients 
One explanation of the reason that available rational-type 
planning models yield significant errors is that in meeting 
the desiderata for generality and simplicity they in general 
do not explicitly include critical driver parameters. Further, 
the functional relationships between the significant drivers 
are not evaluated as regards various flood regimes, e.g. 
surface dominant regions, subsurface dominant regions, etc. 

It would be possible in theory to utilize a more recent genera- 
tion of models, the so-called parametric models, as planning 
models by introducing the data pertaining to the desired re- 
currence interval. This approach presents three problems: 

1. Current parametric models were devised for gaged 

watersheds. They thus require a feedback correction 
through streamgage data of at least a few, perhaps 
five years, duration. 

Current parametric models arO not optimally structured, 
for good and valid reasons, to exploit to the fullest 
the capabilities of remote sensing. 


2 . 
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3. They are quit e-- expensive irf terms of computer time 
and require-reiatively large computer facilities, 

not„ available' to most of the “grass roots" users. 

.1 

It.: -is these drawbacks which th^^ present effort is intended to 

alleviate or obviate. Specifically, this effort is directed 

toward determining which remote sensing observables are most 

. . * • .- * • “ 

important to the planning of watersheds vis-a-vis peak flow. 
Concomitantly, it seeks to Identify the regions in which 

• , . * * ‘ . V /. 

surface parameters .dominate the hydrologic processes, in- 
order to test the remote sensing ' model hypotheses. 
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2.0 GENERAL APPROACH 


With’: the advent of remote sensing, the hydrologist has avail- 
able a practical tool for developing a new data series which 
can be input into planning models. In principle, it is now 
possible to receive synoptic and repetitive information about 
the watershed and to develop therefrom a sequential profile 
of the vegetative cover, impervious area, potential infiltra- 
tion and soil humidity from past potential evapo transpiration. 
■■•The new data eleraent§,are potential drivers for planning models 
their incorporation can be expected to reduce the variance of 
the error. 

The thrust of this effort was -concentrated in four areas: 

1. Determination of driver phenomena, 

2. Construction o;ra' generalized. hydrologic planning 
model, primarily involving "observables..*' 

X 

3. 'Verification of the model with data from existing 
watersheds , 

4. Identification of the role of remote sensing. 

The relationship of each of these tasks to the others arid 

\ I 

their roles in the overall effort are shown in Figure 7- 

<\ 

The ;first''' step in the approach was to attempt to isolate those 
phenomena which "drive" hydrologic plaiining models. A driver 
for a planning model for., peak flow is defined as any watershed 




HYDROljOetC PROCESSES 



nrom OauQiKl Raeottfi 

HCOSYSTEiViS 
IMTERNATIONAL INC. 











- 1 < 9 » • 


condition or parameter which^ when modified, causes a signif- 
icant alteration in the peak runoff. ' 

The driver concept was used to validate three hypotheses: 

1, Not all the physical characteristics of watersheds 
contribute significantly to runoff rate or volume.. 
Therefore, the drivers of peak flow events constitute 
a relatively . small group, of watershed parameters. 

2, Drivers may be spatially and/or temporally variable. 
Much of the uncertainty of current hydrologic planning 

t 

models is due to their inflexibility in accommodating 
natural, and. man -induced variations in watershed 
conditions . 

3, “•’A'fter ascertaining the extent, of influence of each 

driver, it is possible to neglect those phenomena 
'whose contributions to the desired end result (peak 
flow) are minor, i.e., cause an error smaller than a 
preassigned proportion of the total peak flow. 

Once the important drivers were identified, their interrelation- 
ships were investigated. This portion of the effort proceeded 
along two lines. First, the mathematical relationships among 
processes and drivers were explicitly determined, in order to 
develop an analytic model. Next, analog computer equivalents 
were developed. The analytic model was primarily applied to 
the examination of runoff "and’ Its correlation to areal rain- ' 
fall, while the analog representation was used to test the 
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sensitivity of hydrologic processes to changes in various 
basin parameters. 

The proposed verification of model results required the 
assemblage • of a data base- of physical characteristics and 
rainfall/ rurib'ff records extending' ovdr' a large number of 
existing watersheds. ' The advantages -of having' access to 
this -data' base were twofold: . 

1, ' Real information was available during the formulation 
: *v of-the -raodei, • to test the -accuracy of its results. 

;2‘. Since', the - drivers' are"- spatially variable', a large 
ic..:. geographically dispersed sample of basins' facilitated 
• ; c''.'the -determihatlon of which physical ' conditions, and 
therefore which hydrologic processes, were dominant 
in each- -location. - 

Ultimately, the results of this effort should be applicable 

!* r: ‘ I ^ 1*1 , . . * * ' O * * ' • • . ' V 

to the construction of a ’’modular" model, capable of being 
tailored easily to the area under study . . Referring again to 
Figure 7* the A, B, and C versions of the model represent what 
‘are expected tc become --unique combinations, ■ each -including 
• the minimum' numb er of modules to ensure accurate resui'ts. 

■ • •. > 'T * ■ -S-'-J ' ' ^ • ' •'■ - ■ • 

To assemble the required data base, a geographically repre- ^ 
sentative sample of basins was taken from a .group of experi- 
mental watersheds operated by the Agricul'tural Research Service 
(ARS) for the purpose of studying the effects of agricultural 
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practices upon hydrologic processes. 158 watersheds were 
selected from 30 locations in the United States, as shown 
in Figure 8. The sample included substantially all ARS 
'basins with area greater than 100 acres (40.5 hectares),' 
since benefits from improved peak flow prediction are' small 
in watersheds of lesser area. 

Finally, the effort involved the description of the role of 
remote sensing as it relates to hydrologic modeling. Particu- 
lar emphasis was placed upon identification of the areas where 
ERTS-type data is of .maximum benefit and specif ica'tion of the 
method by which specific parameters may be remotely sensed. 
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3.0 INVESTIGATION OF DRIVER PHENOMENA 

/ * 

Figure 9 illustrates the mechanism by whic^i rainfall becomes 
■peak runoff. Th-e figure depicts a watershed of area A, subject 
to rain at constant rate, the runoff from which is measured at 
the watershed outlet^ If th^.- watershed were impermeable and 
the rainfall duration d-ufficiently Ipng^, the outflow would be- 
come constahfat the instant when runoff rate began to equal 
..rainfall rate. Therefore, the peak outflovf rate would equal 
the product of the rain rate times the watershed area. The 
time elapsed from the moment the rain starts to the moment at 
which this peak occurs is termed the time of concentration of 
the watershed. 

A ■' 

Assume, however, that the watershed allows a portion of the 
water to ‘be abstracted (by infiltration, interception, evapora- 

i ‘ 

tion, etc.). Then the peak outflov; rate will be lower; and may 
be determined byr 


^eak ” 

Where : 

P a rainfall rate 

• 

I = abstraction rate 
A, =r .watershed area 
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In practice, as we shall see, I and P are not as simply 
related as the expression above implies. All that can be 
said at this point is that the peak runoff will increase as 
rain increases, and decrease as the abstraction increases, 
but not necessarily in direct proportion to either factor. 

Figure 10 graphically illustrates the principal interrelation- 
ships among hydrologic processes. Table 2 supplies a brief 
definition of each process, while Table 3 lists the potentially- 
important drivers for each process. Note that, “potentially 
important" is only a qualitative indicator. The purpose of 
the discussion which follovxs is to convert this qualitative 
into a quantitative measure of each driver's importance. 

3'.1 ' Precipitation 

Since precipitation is the source of the direct runoff, its 

t 

properties are the basic input to any hydrologic model. 
Hydrologic planning models are concerned with the statistical 
distribution of peak events which are likely to recur within 
a specified number of .years (the 25-year flood, the 50-year 
event, etc.). Therefore, we are concerned here with the 
statistical properties of the rainfall, -e.g. the rainfall 
rate-duration relationship. 

The statistical properties of rainfall, of interest to plan- 
ning models, are: 
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SOMMARy DESCREPnONS OP Hm)LOGIC PK3CESSES 


ffiDHDLOGIC PROCESS 


DESCRIPTION 


Interception 

Kfcdsture caught and stored on plant leaves and 
steins or other inpermeable objects; eventually 
evaporated back into the atmosphere. 

Infiltration 

Downward roovsnent of water from the surface into 
the soil. 

A) Interflow 

Lateral subsurface water movement toward stream 
channels. 

B) Percolation 

Downward movement of water throu^ soil to 
groiindwater (area where pores of soil or rock 
are filled .with water) , 

C) Base Runoff 

' Water from interflow and. percolation which moves ' 
undergroxond' to the channel. 

’ Ev^otranspiration 

T^ard movonent of xvater in gaseous state fron 
the surface. 

A) Evaporation 

B) Transpiration 

Movement of water through plants to the. 
atmosphere. 

Precipitation Excess 

Retention of excess rainfall in surface depres-. 
slons. 

A) Depression 
Storage 

B) Surface Plow 

Uninfiltrated \\’ater which flows over land 
surface to stream channels. 

C) Channel Flow 

Plow of water in mtural channels'. 

Qbtal Runoff 

Sum of runoff from underground processes (base 
runoff) and overland flovr (direct runoff) . 
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TABLE 3 

PCOENTIAILY IMPOFeTANT DRIVERS AS MLAIED TO HYDROLOGIC PROCESSES 


hydrologic process 

PRINCIPAL DRIVERS 

Overland Plow 

*Slope 

*Eoughness of Soil 
& Cover 

^Drainage Density & 
Pattern 

Infiltration 

*Soll Permeability 
^Antecedent Soil 
Moisture 

Soil Moisture Capa- 
city 

A) Interflow 

*Soil Permeability 
Subsurface Moisture 
Gradient 

• *PlOTf Length, Slope 

B) Percolation 

*Soil Permeability 
Subsurface Itolsture 
Gradient 
Soil Depth 

Evapotranspiratlon 

Tenperature 
*Antecedent Soil 
Moisture 

*Soil Permeability 

A) Evaporation 

B) Transpiration 

Tenperature 
* Solar Radiation 
*Vegetatlve Cover 
*Antecedait Soil 
Moisture 

Depression Storage & 
Detention 

*Depression Density 
*Cover Retention 

Interception 

Duration of Rainfall 
Intensity of Rain- ’ 
fall 

*Cover Conpositlon, 
Age, Density 


SECONDARY DRIVERS 



^Vegetative Cover 
*Slope 

Water Turbidity 
Tenperature 



Water Turbidity 
Wind 



*Slope . 



* indicates factors Which are potentially remote sensing observable. 
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1. The magnitude of the peak rain event varies with 
geographic location, as shown in Figure 11. For 
the United States, the range of variation is approx- 
imately 3.5:1, if unusually dry areas are excluded, 
such as Death Valley, for example. 

2. The rate of the peak rainfall event varies as a 
function of the length of the event, as shown in 
Figure 12. The longer the event, the lesser the 
rate. However, the mass, i.e, the total amount of 
water precipitated, grows with increasing event 
duration. 

3. For a given region, the rain recurrence-durat ion- 
intensity relationships follow the empirical 
relationship: 



Where : 

i = rain rate, cm/hr. 

T = recurrence internal, years 
t = rain duration, hrs. 

tti, a 2 

« 3 , d =* constants for specific location 

The constants can be determined from existing rain- 
gage records. .Figure 13 illustrates a typical re- 
lationship, related to the curves of Figure 12. 
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rainfall RATE/ MASS RELATIONSHIP 

DANVILLE, VERMONT WATERSHEDS 
50 YEAR RECURRENCE RAIN 

MAXIMUM RATE AND MASS VERSUS DURATION 
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DANVILLE, VERMONT EXPERIMENTAL WATERSHEDS 
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General Formula 
where, 


<L- 


i = rain rate, cm/hr. 

T= recurrence interval, years 
t= rain durotion, hours 

* constants for specific location 

3.30 

For Danville ; 
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The intensity of rainfall is not uniform within 
the a.rea rained upon. The larger the area,- the 
greater the reduction in rainfall rate as one 
■progresses from -the point of highest rate towards 
the edges of the area. 

Therefore, when considering v/atersheds of relatively 
large extent, a correction factor for rainfall should 
be employed, as indicated in Figure 14. Note that 
the areal distribution of rain is only imperfectly 
known, because the experimental data colle'cted up 
to now on this phenomenon have been scant. Thus 
Figure l4 should be taken as indicative only. 

The time scale of interest in the abscissa of the curves of 
Figures 12 and 13 can be approximated by computing the con- 
centration. time from existing formulations. The result is 
only approximate, but sufficient to yield a gross calibration. 
The computation, whose details are reported in Section 5, is 
summarized in Figure 15, from v/hich it can be deduced that the 
times of interest for the watershed sizes, shapes and slopes 
contained within the 158 ARS test watersheds range from perhaps 
10 minutes to upwards of 5 hours. 
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3*2 Subsurface Processes 

3.2.1 Subsurface Abstraction 

A fraction of the rainfall reaching ground level will penetrate 
or infiltrate into the soil. This fraction — or at least, the 
portion thereof that does not evaporate — eventually becomes 
runoff via the subsurface flow mechanism. However, as will be 
shown, the runoff from the subsurface component is almost always 
considerably delayed with respect to the runoff from the surface 
component. Thus, except for very special cases, the portion of 
rainfall which is abstracted to the subsurface does not con- 
tribute to peak flow events. The fundamental role of the ab- 
straction is to reduce the magnitude of the peak event. The 
crucial question, as regards hydrologic planning models, is 
thus: How- much of the peak rainfall is abstracted? 

There are three major drivers of subsurface abstraction: 

° Soil permeability 

® Soil storage capacity, a combination of porosity and 
soil depth • ' 

° Antecedent soil moisture 

Permeability is conditioned by the degree of resistance to 
fluid flow through the soil. It is essentially determined by 
soil type and humidity. If is defined as a rate, in terms of 
centimeters per hour of water absorbed by the soil. For a 
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given, type of soil', permeability assumes the largest value 
when the soil is dry; this is known as- the initial permea- 
bility. • As the soil becomes wet, permeability decreases, 
until it assumes a steady state value, .known as the final 
permeability. Figures l6 a and b give initial and final 
permeability or infiltra.tion rates, for the range of soils of 
interest to most hydrologic planning models. 

Porosity, the percent of total soil volume available for water 
storage, does not per se affect the abstraction rate — also- 
known as infiltration rate - but rather affects the amount of 
infiltrated water which may be stored. For a given soil .por- 
osity, the shallower the soil layer above the impermeable layer 
the lesser the -capability of the soil to store water. When 
this layer becomes saturated, infiltration into the soil ceases 
If there are several layers of different porosity and depth, 
the one'.wi-th the least permeability controls the process. If 
the soil layer is very deep with respect to the total rain- 
fall-, it can be essentially considered as an infinite layer. 

In this case, the subsurface abstraction is controlled only 
by the soil’s permeability. Figure l6 c shows that highly 
permeable soils, such as sands and sandy loams, do not retain 
water V7ell because they have low porosity. Less permeable 
soils such as clays, however, have higher rates of water re- 
tention, due to high percentages of porosity. 
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Available moisture capacity is the total amount of water which 
a given depth of a given type of soil can store. It is a 
measure of the maximum- amount of water v/hich can fill all the 
available pore space. Figure 16 d lists the available moisture 
capacities for various types of soils. Soils with low infil- 
tration rates and high porosities have larger values of avail- 
able moisture capacity than do soils with high infiltration 
rates and loiir porosity. 

Antecedent moisture is the amount of water already present in 
the soil at the beginning of a rain event, and has the effect 

I 

of reducing available moisture capacity and infiltration rates. 
Its net effect upon the peak runoff event is to increase. its 
magnitude, the more so the greater the antecedent moisture. 

What this means, in effect, is that a peak rain event of a 
given recurrence does not necessarily generate a peak runoff 
event of the same recurrence: the 50-year flood does not 

necessarily correspond to the 50-year rain. It can, in fact, 
be generated by two lesser rains - for example, 20-year rains - 
occurring sufficiently close together in time. 

In- hydrologic management models, antecedent soil moisture is 
handled- by •computing 'the. net difference .between the- water in- 
filtrated into the soil and the water lost through- evapotrans- 
piration, over an interval of time preceding the event under 
study, which Interval can range upwards of a few months. 
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In.,hydrologic planning models, the antecedent soil moisture 
must be evaluated by using statistical methods yielding the 
probability, in. a given region, of having a_ certain sequence 
of peak rain events, separated by specified time intervals. 
This investigation is reserved for future phases of this re- 
search effort. 

Factors such as slope, water turbidity, and temperature also 
have an effect on subsurface abstraction, but the magnitude 
of their Influence is negligible when compared, to the influ-, 
ence of permeability, porosity and antecedent moisture;. 

Figure 17 shows the effect of soil permeability upon infil- 
tration for two ’’infinite” soils in Prince George’s County, 
Maryland,’ derived' from an an^og computer slraulatioh. For 
a given rain e’vent, the low' permeability of one of the soils 
(Lednardtown silt loam) causes reduced infiltration and in- 
creased runoff mass; the high permeability of the other soil 
(Groom gravel loam) allows k higher rate’ of ' infiltration, and 
a consequently lower runoff mass. Therefore, a watershed com- 
posed of deep Lednardtown silt loam would discharge more run-' 
off and have a shorter time of . concentration (i.e.., would pro- 
duce peak runoff faster) than v/ould a similar watershed made 
up. q.f .Groom gravel loam. 

Figure 18 illustrates the importance of soil storage ’ capacity 
(porosity and depth combination) to infiltration and runoff. 
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FIGURE 18 INFLUENCE OF SOIL DEPTH UPON RUNOFF 
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Here, the finite depth of the upper layer of soil limits 
total available storage. When the upper layer becomes 
saturated, there is a sudden drop in the effective infil- 
tration (or abstraction) rate because there is no more stor- 
age available. The water cannot sink as readily into the 
soilj consequently, the runoff increases. 

, -f 

Figure 19 demonstrates the effect of antecedent soil moisture 
on infiltration and runoff. Initially dry soil accepts water 

f 

readily; therefore, time to beginning of runoff is high and 
runoff mass is low. Hov/ever, when there is a large percentage 
of moisture already in the soil, less water inflltrat.es; the 
runoff begins earlier in time and mass is larger. 

Although a great deal of research has been conducted in the 
area of subsurface abstraction, no consensus presently exists 
as to which mathematical formulation best describes the- process. 
Table 4 lists the most widely applied formulations for the 
effective infiltration. These include the effects of perm- 
eability, storage capacity and antecedent soil moisture. Most 
of the equations basically describe an exponentially decaying 
function which finally declines to a finite value equal to 
the saturated conductivity of the soil. The Holtan equation, 
graphed in Figure 20 using actual data from the Danville, 
Vermont watershed, is illustrative. 
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where K,a « constants, empirically determined; no physical significance 
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where L ® depth to wetting front 

f L = matric potential of wetting front 
K =* soil permeability 
f a soil pcrosity 

b)fp® Kg 1 + (M^ * S/F) 
i « If + 

where fp® Infiltration capacity 

Kg® saturoted cr^uotlvity If 

M<j® initial moisture deficit, vol./vol. 

S ®caf8llQry suction. Ins. 

F ® cumulative infiltration 

Phimp 
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Q* At * Bt 
Qal/eAr*^^+B=i 

where A, B® empirically determined constants- no physical relation 
Holton 
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where SraQvolloble porosity® So*Fp 
k ® v^tative cover constant 
fc« final infiltration rate 
Pp® cumulative infiltration to time of fc 

numerical constants are empirically determined 
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Principal Infiltration Formulations: 


Hotton; Overton ( Huggins a Monk Modification) 
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where A, B “constants characteristic, of a given soil a antecedent condition 

S ® storage potential of soil within infiltration control zone® 
total pcs'osity-cntec^iQrjt soil moisture 

Tp “total soil porosity 


Richards 
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In spite of differences in the mathematical expressions, 
computations for two of the most widely used formulas 
(Holtan and Horton),, shown In Figure 21, with the same 
initial and final conditions, indicate that the differences 
in the. results are relatively small. A lesser-used for- 
mulation, that of Green and Ampt, yields a different, end 
infiltration when normalized to the same initial conditions 
as the other two formulations. Conversely, the Green and 
Aupt formula yields different initial conditions when nor- 
malized to the same end-conditions . 

\ 

Approximately • similar overall results can be obtained by 
deciding, case by case, whether it is more representative 
to favor initial or terminal conditions. For short rains, 
for example, initial conditions control. 

The Holtan formulation was tentatively and temporarily 
selected for use in the model on the basis of its suit- 
ability tp remote sensing: 

I = a. GI .(Sa - if 

V/here: 

I = infiltration rate 
a = average vegetative cover factor 
GI = growth index 

Sa = average available soil moisture capacity 
I = cumulative subsurface abstraction 
If = final infiltration rate equal to saturated 


conductlvitv 
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There are two advantages in using Holtan's equation: 

1. It includes several of the physical factors which 

affect infiltration: surface typq. .(determined by 

j*'»‘ 4?^ .-.'K • ■••• • '' > 

■ ■ vegetative cover), permeability (included in the 

If term), and available soil moisture capacity. 

2. All of the factors included can be determined or 
at least inferred from surface observations, with 
a and GI being directly measurable. Values of a. • 
can be determined by Identification of the existing 
cover and GI values in a range of 0 - 1.0 represent 
the maturity of the agricultural cover or crops 
present, with a value of 1.0 equalling maturity. 
Storage capacity and final infiltration rate are 
both Inferrable from soil type and soil, moist.ure. 

In conclusion, infiltration rates, given earlier in Fig- 
Tore 16, vary from as low as 0.13 c'in/hr to as high as 25 
cm/hr, rain rates vary within the same range. Thus ef- 
fective infiltration rates are a svA)stantial percentage 
of total rainfall rates. Infiltration is, therefore, a 
critical driver in any hydrologic planning model. 

3.2.2 Percolation and Interflow 

The remaining subsurface processes are percolation and 
Interflow. Percolation is the process by which infiltrated 
water drains deep into the soil to eventually become ground- 
water, Interflow is the process by which infiltrated 
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water flows laterally, eventually becoming surface runoff. 
To have significant impact upon hydrologic planning mod- 
els, percolation and interflow must be of sufficient mag- 
nitude to significantly alter the runoff mass. Also, the 
water Involved must reach the watershed outflow point in 

time to contribute to the runoff peak. 

1 

In '-real watersheds, with possibly rare exceptions, neither 
of ^the above requirements is met. Rates of percolation 
and 'interflow are lov7 to negligible compared to those of 
overland flow, as computed in Table 5 ‘and summarized in 

\ S 

Figure 22. This difference in speeds, is due partly to 
the fact that during a peak flow event, the pressure of ‘ 
water flowing in the channel is opposed to the flow of 
subsurface water, thus retarding the Inclusion of sub- 
surface water in the outflow volume. 

3.3 --^Evapotranspiration 

Evapotranspiration combines two processes: 

1 ‘ , ' 

Evaporation - the loss of water to the atmosphere. 

\ ^ * 

„ w 

• ^ ' 

° Transpiration - the process by which water is drawn 
out of the soil and transpired by vegetation. 

Evapotranspiration can account for direct Tosses of pre- 
cipitation and for losses of stored, infiltrated, and in- 
tercepted water. 
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Thei*® are five major drivers of evapo transpiration: 
Temperature 
*0.. Solar ‘Radiation 
o Antecedent -Soil- Moisture 

I « 

® Soil Permeability 

- : 

i -0 Vegetative Cover 

Temperature and solar radiation are related since they 

'?• 

both function in determining amount of heat present in the 

air. When the air temperature is high, evapotranspiration 

■ ^ 

rate is also high since available moisture evaporates 

/ 

quickly. When the air is cool water evaporates more slowly 
and evapotranspiration rat'^ is lovr. 

Evapo t inspiration is also affected by antecedent soil 
moisture, since the process can continue only long as 
moisture is present in thevboil or on the surface. Once 
this moisture is depleted, evapotranspiration must nec- 
essarily cease,' even if other drivers are present. 

Soil permeability determines '-.the rate.'at^--Which water can 

move upward through the soil, just as "it influences the 

• • ■ . 1 / ^ 

rate of infiltration. Soil permeability, which is deter- 

^ t 

mined by soil type, and rate of evapotranspiration are less 

for clay soils which resist movement of water through the 

soil than for sandy soils which have low resistance to fluid 

movement. •• 

ORIGINAL PAGE IS 
OP POOR QUALH^^ 
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Vegetative cover is Important because plants transpire 
differently. In an area where there are not plants, evap- 
otranspiratlon involves evaporation only and its rate will 
be determined solely by the rate of evaporation. However, 
in an area where vegetative cover is iieavy, rate of trans- 
piration must also be considered. 

Other factors also influence evapotranspiratio'n. Wind," 
for example, functions in removing "moist air from above 
vegetation and v/ater turbidity or inuddiness slows down the 
rate' of evapo transpiration. However, the influence of • 
these factors is quite small compared to the effect of’ 
temperature, solar radiation, antecedent soil moisture,- 
soil permeability, and vegetative cover. 

To ascertain the importance of evapotranspiratipn relative 
to other parameters, average rates vjere calculated from 
the ARS .test watershed data, using^the Thornewaite equa- 
tion: 

✓ 

Avg. ET = .022 (l^\^ 

... , TE< 

’ Vihere: 

' ,* * 

Avgl ET = average evapotranspiration 'rate. Cm/hr. 

t'= mean monthly 'temperature, "‘^C ' ' 

= S (t/5)^*^^^ for each month 
a = .000000.675(TE)3 - .0000771(TE)^ + V01792(TE) 

+ .49239 



Results are presented, in. Table 6. The maximum practical 
value is approximately 0.4 cm/day. linien compared to rain- 
fall and Infiltration rates, this indicates that evapo trans- 
piration occurring during peak events is unimportant to hy- 
drologic planning models. 

However, evapotranspiration preceding peak events is impor- 
tant, because its effect is that of depleting antecedent 
moisture, thus increasing the subsurface abstraction. For 
example, in regions where rate of evapotranspiration is 
high and the rainless period of long .duration, antecedent 
moisture will be greatly reduced.- For example, at an av- 
erage rate of 0.4 cm/ day, an area which experiences two 
storms spaced 8 days apart will lose 3.2 cm in upper soil 
layer water content betvreen storms, significantly Increasing 
storage capacity in the soil. Therefore, although evap-' 
otranspiration is minimal in its influence upon runoff, it 
is Important in its effect upon soil moisture content. 

This importance is a function of the available storage cap- 
acity. In shallow soils, for example, or in deeper soils 
with low porosity, the subsurface abstraction is a smal- 
ler-, fraction of the precipitation than for deep, porous 
soils. Hence, for the form.er type of soils, the antecedent 
moisture- plays a relatively minor role. As. we shall see, 
it is exactly these soils which are the most amenable to 
the construction of hydrologic planning models based upon 
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remote sensing techniques. 

3.4 Interception 

Interception is the process by v/hich rainfall is prevented 
from reaching the ground to become runoff or infiltration. 

It is dependent upon duration and intensity of rainfall, 
upon the species, composition, age and density of vegetation, 
.and upon the season of the year and the geography of the 
region. 

Two factors limit the importance of interception to the 
hydrologic modeling of peak fTov; events. First, intercep- 
tion stores only small quantities of water. For example, 

in a basin completely covered by deciduous forest, and 

1 

thus of highest interception capacity, interception v;ill 
claim about 0.35 centimeters of a 3*5 centimeter rain. 

Thus, in practical cases but little of a peak rain will 
be subject to Interception storage. ■ Second, the time per- 
iod required to fill interception storage is negligible. 
During a 1.-8’ centimeter per' hour rain, for example, only 
approximately 15 minutes are required to saturate inter- 
ception capacity. Rain intensities of interest to plan- 
ning' models are generally on the order of three to five 
times greater; saturation times will be proportionately 
smaller . 
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3 « 5- Depression Storage 

Depression storage is similar -to Interception; it is the 
process by which vmter is retained at the surface by reces- 
ses ^ or shallow small pools, in the soil. It is determined 
primarily by depression density and retention of cover. 

Like interception, depression storage is- not - a continuing 
process; it ceases when a fixed capacity is reached. . Typ- 
ical saturation capacities are 0.5, 0.375, and 0.25 cen- 
timeters for sand, loam, and clay soils, respectively. 

In the Washington, D.C. area, for example, only approx- 
imately h% of the 50-year 1-hour rain would be stored by 
depression storage. 

Depression storage is thus of minor importance to the mod- 

f . “ ‘•f • */* 

eling of peak flow events. 

3.6 Overland Flow 

'During ’a' peak event, the rate at which the excess rainfall 

oVer and above the subsurface abstraction contributes to 

streaniflow is determined by the overland flo>r process. The 

process is shown schematically in Figure 23. The excess 
^ ^ % 1 ' 
rainfall, i.e. the rainfall less the subsurface abstraction, 

flows over the surface until it encounters a. natural chan- 
nel. It is then conveyed from the higher-order channels 
into successive lower-order streams. In the channel, the 
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.water ,g.enef*aliy-, speeds .up since it travels faster in the 
channels than over -the surface, it v/ould be expected that 
the time of concentration, i.e. the time required for the 
water to travel from the uppermost reach of the watershed 
to its outlet, would be shorter the more the number of 
channels per unit area of the v^atershed. The total length 
of the channels divided by the area of the watershed is 
known as the drainage density, it has been indicated in 
the literature, and demonstrated in this project, that 
the drainage density plays an important role in deter- 
mining the magnitude of the peak. flow event. 

( ' * 

The rate of overland- flow is governed by the fpllovjing: 

° Watershed slope 
o Watershed cover 
-O Drainage density 

® Channel pattern 

The ground slope determines the magnitude of the grav- 
itational force which induces flow of the water. The type 
and density o.f ground cover determines the frictional force 
opposing flow. Hence the combination of slope -and cover 
is a determinant of flow velocity. The drainage" density 
defines the average distance over which water .must ' travel 
bdfore reaching the stream. • The channel pattern- establishes 
the-" speed with which water is' carried away by the- channb'l 
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network. In combination with surface flow velocity, this 
defines the time of concentration. 

In the literature, overland surface flow is generally con- 
sidered to be sheet flow. The reason alleged is that the 
depth of flow is small, compared to its width. Implying 
that viscous fluid forces predominate over those produced 
by inertia. Further, overland flov; is often considered 
to be uniform; that is, no change in magnitude or direc- 
tion takes place over the flow length. In this project, 
the assumption of sheet or laminar flow was retained as 
a preliminary working hypothesis, to be analyzed further 
for correctness if vrarranted by significant discrepancies 
between results and predictions from the model. 

Overland surface flow can be quantified by assuming the 
land surface to be approximated by a wide channel. Since 
channel sides have minimal effect upon velocity when the 
channel v/idth exceeds ten times the depth of flow, the hy- 
draulic radius for overland surface flow becomes essentially 
the depth of flov/. 

Overland surface flow eventually reaches a channel; thence 
it is discharged into a sequence of -ever-larger streams, 

.as. it proceeds towards the watershed outlet. The cross 
sectional geometry of natural channels is a major deter- 
minant of flow rate .within the channel.- Figure 24 graph- 
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Ically demonstrates its effect upon channel flow. As the 
width to depth ratio increases, the exact channel geometry 
becomes less important: the hydraulic radius can more and 

more be approximated by the depth of flow. 

The particular cross sectional shape that evolves naturally 
in a stream is that which best balances the transmission 
of water and the structural stability of the banks, in 
channels with firm banks,^ for example, the shape most ef- 
ficient to flow is parabolic. However, the cross sections, 

1 

of natural streams are not generally parabolic in shape: 
but -appear to be trapezoidal, at 'least in the stream's 
straight, sections. Empirical evidence has shovm that chan- 
nels- in the dry vrestern states are usually wider and more 
shallow than those in humid areas. 

The shape taken by channels is determined by the peaic ev- 
ents- of a stream rather than by its average flows, pos- 
sibly because during peak events, the forces acting upon 
stream beds and banks are the highest. Typically, U.S. 

rivers carry less than their mean flovj 60-75^ of the time, 

’ ^ \ ♦ \ 
and less than half, of the mean jfflov? about 25^ 'df "the; time . 

Average annual discharge -1*1115 the channel to only approx- 
imately 1/3 of its bankfull depth- Bankfull discharge oc- 
curs about once every 1.5 years; a flood plain is inundated 
to approximately 1.8 times the bankfull depth of the chan- 
nel once in fifty years . 



The channel cross section changes .as a function of lbs 
position along the stream axis. Progressing in the down- 
stream direction, for example,, -discharge increases, causing 
stream geometry to alter. The channel's top vridth in- 
creases approximately as the square root of discharge. Al- 
so, because of higher flows, bed roughness tends to dim- 
inish. The trapezoidal shape becomes progressively more 
rectangular, since stream width increases faster than depth 

t 

; ; 

The Manning equation provides a good description of the 
relationship between the hydraulic parameters, roughne'ss 
and mean channel velocity: 

Velocity (v) = i. R 
m/sec 

Where: 

R = hydraulic radius, m 

S = slope, m/pij ' 

_l/-3 

n = Ifenning’s roughness coefficient, m ' sec 

The hydraulic radius, defined as the cross sectional area 
of flow divided by the wetted perimeter, varies with chan- 
nel; geometry, as shown in Table 7. In natural channels, 
the variations are usually small. Since the top width is 
typically of the order of, or greater than^ 10 times the 
depth, the hydraulic radius for all common cross sections 
mas\-be approximated by the depth of flow without signif- 
icant loss of accuracy. 



Section 


Surface Flow 


Approximoto 


Areo 

A 


Wetted' perimeter 




Hydraulic radius 
R 


xy 


X4.2y 



2x^y * 


3X2+ 8y2 



XAB4^ PmmE'P^-^S C^^N^CRgSi- 




























- 67 “ 


Thus,' Manning’s equation, for both overland surface flow 
and wide channels, becomes;- 

Velocity (V) = .i 

n 

Where; 

d = depth of flo.w, m 
S = slope, m/m 

■ -!/■=? 
n = Manning’s roughness coefficient, m sec. 

The corresponding equation for the flow is; 

. -Plow (Q) = Ay- 

■3 

m-’/sec 

Where: 

2 

A = cross sectional area of flow, m . 

Manning’s roughness coefficient (n) is a measure of the 
frictional resistance of the surface to flow, due-to the 
watershed’s surficial characteristics. The roughness coef- 
ficient is a combination of the effects of various surface 
factors. The value- of n for channel flow is -influenced by 
the following factors: 

;• o Roughness -. • 

Presence .of vegetation 

° Irregularities 

° Channel alignment 


Silting and scouring 



;v:' P< Ot»'st'ruct Ions' 

° Staee and discharge 
Seasonal change 

® Suspended material and bed load 

Similar effects are present in overland flow, although the 
magnitude and importance of each differs. Surface roughnes 
vegetation and seasonal change, for example, are much more 
important to overland=' flow than they are to' channel floxv, 
while silting and scouring, stage and discharge, suspended 
material and bed load are more important to channel flow. 

Manning's n then can be written ,as: 

c . . 

n = mS(n^) 

n = the effective Manning's n 

ri ' 

' = contribution of roughness factors to n 
m = meander factor 

Typical. value;S of: the/ and m are given by Vente Chow, 
Handbook of Hydrology , and are reproduced inTahle-S. 

For example, n for a gradually varying coarse ‘ gravel ' chan- 
nel with minor irregularities, negligible obstructions, 
high vegetation, and minor- meandering would be: 

n = (0.028 + .005 .000.+ .037'^ 1 ‘ 



Typically, n can range from 0.03 to 0 . 29 . 
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Channel cc^itic^s Voluss 

Moterial 

involv^ 

..Eorth • 


X5J52Q 

Rock cut 

0.025 

Fine grovo! 

0-024 

Goarsa groval 

o;o23 

Degree of 
Irregularity 

Ssn^th 


0^0 



fivjsigrate 

0.010 


0.0^ 

Variations of 
Channel cross 
section 

'Gradual 

'^Z 

O.OCKD ' 

Alt^otsf^ ©5co3ien!^ly 

o.oos 

. Alfematisg treqt^y 

0.010^0.015 

Relative i 

Eff^t of 
Obstructions . 


”3. 

ii 

s:soD 

f^inor 

0.010-0.013 

Apprsoiabie 

o.oso-o.cso 

Severs' 

0.040- O.C^ 

Vegetation . 

Low 


0.003-0-010 


aoi 0-0.023 = 

High 

0.025-0.050 

Vsry high 

0.050-0. !00 

Degree, . of 
Meandering 

Minor 

nu 

o 

I.PPO 

Ap^ociot^e 

U50 

•Seysrs 

1.300 
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Typlcal values of n for different overland surfaces are 
given in Table 9* is clear from these figures that 

land cover is significant in determining runoff. Differ- 
ences in surface roughness can cause overland flovj vel- 
ocity to vary up to approximately eight times its min- 
imum value. 

Routing the runoff mass to the basin outlet involves con- 

/ 

sideration of both overland and channel flow. The dif- 
ference between the two types of flow determines the amount 
of water stored on the surface or built up on the chan- 
nel, and, therefore, the depth of flow and its rate. 

3.7 Important Processes and Drivers for Hydrologic Planning 
Models 

Figure 25 summarizes graphically the relative magnitudes 
of the principal hydrologic processes. The ranges of val- 
ues shown in Figure 25 represent the ranges encountered 
in the 158 ARS test watersheds. 

It is clear that the primary hydrologic processes which 
dominate peak flow events, and whose representation, there- 
fore, must be included in a generalized watershed planning 
model, are rainfall, infiltration or subsurface abstrac- 
tion, and overland flow. Secondarily, according to the 
region under consideration, the statistical behavior of 
antecedent moisture condition as it relates to peak rain- 
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Tobls 9 

mmim's roushmess coefficient fc^ overland 

FLOW FOR VARidB SJRFACE TYPES 


Watershed Surfoce 

Manning’s “N" 

Smooth Asphalt 

0.013 

Concrete (Trowel Finish) 

0.013 

Rouah Asphoit 

0.016 

Concrete (Unfinished) 

0.017 

Smooth Earth ( Bare) 

0.018 

Firm Gravel 

0,020 

Cemented Rubble Ma^ry 

. 0.025 

Pasture (Short Grass) 

0.030 

Pasture (High Grass) 

0.035 

Cultivated Area (Row Crops) 

0.035 

Cultivated Area (Field Crops) 

0.040 

Scattered Brush, Heavy Weeds 

0045 

Light Brush and Trees (Winter) 

0.050 

Light Bnsh and Trees (Summer) 

0.060 

Dense Brush ( Winter) 

0.070 

Dense Brv&h ( Summ^) 

0.100 

Heavy Timber 

0.100 

idle Land 

0030 

(Brass Land 

0.032 


ECOSY. 

INTERNATlOt 



Rate, cm/hr. 




50 Yr. rain- range over 150 Experln^ntai 

watersheds 

I Infiltration- '• » « »• h 





I Evapotranspi ration m 

AJJUiXULAXjLJLJLLA.^^^^ ... ... C0^^ 


Interception 

Depression Storage 


12 16 
Time, hrs. 


■fall events must- be Included. 

. . ^ ^ . 

Rainfall is the principal causative factor defining the 
magnitude of the peak flow. Its -important characteris- 
tics are the recurrence statistics, determined by empir- 
ical correlation of regional rainfall records. 

*.!|fnfiltration -governs the portion .of the ra.infall which 
...contributes to the direct runoff- peak-. It-.can.;be evaluated 
; from . wat.ershed soil records^ abundantly available. 

The overland flow process and channel flow determines 
the timing of the peak. , The timing in turn, determines 
thet .rainfall, rate an,d mass for a_^ given- recurrence fre- 
quency, .and .hence determines -the peak ;flow. The overland 
flow can be modeled from knowledge of the surface charac- 
teristics of the watershed, which are directly amenable 
to remote sensing. 

•The key drivers of peak flow, in addition to ’rainfall sta- 
• tistics ..are: . .• -.r 

1. Soil Permeability - high permeabilitie,s -mean high 
acceptance of water and smaller runoff mass, 

.2. • Soil Water Capacity - a soil having a greater water 
water capacity will retain -more rainfall and produce 
less' runo-ff . 

3 .' Antecedent Soil Moisture - as soil moisture rises, 
the soil becomes saturated, slowing infiltration 
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rates, reducing total soil moisture capacity,- -and - 
increasing the runoff volume. 

4. Slope - flow velocity varies directly but non-lin- 
early with slope. 

5. Surface Friction - velocity varies inversely (and 

,: r • j • vjI 

non-linearly ) with surface friction. 

= 6 Drainage Density • arid Pattern - -^defines the relative 
.fA?'-- -tiistances that water will flow overl'and" and ‘in the 

dhannel;' in combination with slope aft’c^'^surface fric- 
tion, defines concentration time. 

A first-cut quantification of- the sensitivity of the run- 
off volume, or of ■= -hydrologic quantities -impacting runoff 
'volume, to variations In 't'hese' drivers -is- -presented in 
Table 10;. ' • ‘ 

Areas which are dominated by surface processes (i.e., pro- 
duce the most surface water per volume of precipitation) 
-will iri general derive the -most benefits from hydrologic 
planning models. Also, regions which are surface dominated 
are the-'-best- served by remote sens'irig-. 

The influence of temporal variations of watershed par- 
ameters - for example, caused by urbanization - well amen- 
able to remote sensing, is more prominent in surface dom- 
inated regions. 
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Table* 10 Sensitivity of Runoff to the 
Principal Drivers 


Driver 

Practical Range 

Computed 

Effect Upon Runoff . 

Slope 

0.01 to 1.0 m/m 

lb k^increase in 
flow rate ;( M^onning) 

Surfoce Friction 

0.01 to .lbm“*^*sec 

- . . • • r- 

10 X decrease in 
flow rate (Manning) 

Drainage Density 

/ 

i/loo to i/iOdoo 

m/m* ' ^ 

Up to approx (Ox 

>:d 0 crease in flow ' 
velocity 

Antecedent So H ' 
Moisture . 

«•/ . ' ^ 

1 

0 to 100% of ■ ' 

. available capacity 

Up to' 2x dsjrease 
;ln runoff volume 

I 

SoilWater'^ 

Capacity 

K/f-. — 'f- .. 

‘ ’ 0 to 60% of ' ■ ‘ 
upper soil layer 
volume 

Up to 2x decrease" 
in runoff volume 

Soil Permeability .. 

0.Q02.5j to .25 m/m; 

up to iOx decrease 

in rate of pr^iplta- 
tlon excess butit^p 

't * * * - 
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Data' acquired from the 158 test watersheds, a summary of 
which is contained in the Appendix, provides the basis 
for partitioning the United States into areas possessing 
either surface-dominated or subsurface-dominated hydro- 
logic regimes. From rainfall and runoff data, average 
annual precipitation (P) and average annual discharge (Q) 
were calculated for all test basins over the period of 
record. Q divided by P yields a first-cut measure of the 
propensity of a viatershed to discharge. Figure 26 gives 
Q/P numbers for the regions in which the test watersheds 
are located. On the basis of these calculations, an in- 
itial partition of the United States into 3 categories 
of hydrologic regimes appears as in Figure 27. The three 
regions are: 

1. Heavily Surface Dominant - Where the percentage 'of 
rainfall to runoff significantly exceeds the per- 
centage of rainfall to infiltration. 

2. - Surface Dominant - Where more rainfall runs off 

than infiltrates. 

3. Subsurface Dominant - V/here more rainfall Infiltrates 
than runs off. 

It is interesting to note that the regions which are sur- 
face dominated are also those which have historically ex- 
perienced the greatest flood damage. 



AGRICULTURAL RESEARCH BASINS loo- 100.000 acres 


Average Annual Runoff t Average Annual Precipitation 



. 100 - 1000 ACRES 
o 1000 - 10,000 ACRES 
* 10,000-100,000 ACRES 




INITIAL RARTITION 

BASED UPON 

AGRICULTURAL RESE^GH GAGINS 100-100,000 ACRES 



I 

00 

1 


Surface Dominant 
S3 Subsurface Dominant 
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^.0 GENERALIZED HYDROIjOC-IC PLAIMENG M3DEL 

This project attained two milestones tovrards the achievement of a 
hydrologic planning model: 

1. An overall framevrork was formulated. 

2. A number of segments, or modules, were constructed, each 
modeling a distinct hydrologic process. 

V/hat remains to be done is to connect and Integrate the various 
modules into a single model. 

The modules consist of, and the Integrated model is eJ^ected to 
consist of, analytic expressions suitable for digital or analog 
catputer progranraing. 

The analytic procedure was supplemented by an analog cotputer simu~ 
lation, which was oriented primarily at estab lishin g the sensitivities 
of the runoff to variations in the inportant watershed parameters: 
soil permeability, antecedent soil moisti^e and total storage . capa- 
city (soil depth). The advantages of the analog sinulation were 
found to be: 

1. Flexibility and ease of variation of rates and magnitu<^s 
of the physical parameters. 

2. Adequacy in representing the physical phenomena, yet with 
relatively small conputer hardware. 


Figure -28 is a block diagram of the analog conputer model developed. 




figure 28 ANALOG CIRCUIT FOR 

CONSTANT RAIN WITH EXPONENTIALLY DECAYING INFILTRATION ; CONSTANT ET 
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It lends itself to simulating rain events of varying characteristics: 
rates, durations, and rate variations vjithin the rainfall, period. It 
also simulates the subsurface abstraction process \d.th varying parar 
'meters: ’ initial and final infiltration rates, infiltration decline 
rate. Derived parameters v^ich can be simulated are soil depl^, soil 
storage capacity, antecedent moisture. The output is the.tot^ 
ijie overland flow conputation is performed analytically, althou^ ^ 
analog simulation is quite feasible and is contenplated for future 
phases' of this effort. 

Sinc^, as has been shown, the effects of interflow, percolation,' 'arid- 
depression storage are minimal during peak-, events., = they are not ih^'’- 
cluded explicitly. They can, however, be factored in through potentio- 
meter P8 of Figure 28, as a constant small rate. Abstractions from 
evapotrarisplration can be factored in similarly. Although, as tes 
been in the preceding section, the abstraction friom evapotrans- 

plration is 'generally snail mth respect to rainfall and to subsurface 
abstractions, there may exist combinations of meteorological and . 
physical watershed parameters vihich may warrant its inclusion. For 
example, evapotransplratlon consideration should be included when one 
of the following conditions exist: the watershed to be modeled is 

contained in an area of high evapotransplratlon potential (dry sunny 
climate.); it contains highly evaporative vegetation; the peak recurring 
i^ns are riot .very intense; -the times of concentration are' long' - 
(large watersheds v/ith high surface frdotion) ; or- the: precision desired 
in .the model’s output is very high. In most practical cases, however. 



.,-82- 

:^;he ewapotranspiration component can be neglected. 

Consideration of evapotranspiration becones iiiroortant when the peak 
flow event is produced by multiple sequential rainfall events. In this 
case, if the statistical distribution of the inter-event intervals 
shews that they can be relatively long, evapotranspiration becomes 
lnportant in determining the statistics of antecedent humidity. Al- 
though multi-storm analysis was not performed in this phase of the 
project, analog simulation circuits were devised to provide the 
capability for such a study at a later time, llie circuit is depicted 
in- figure 29. It simulates a sequence of rain events of arbitrary 
intensity, duration, and inter-event period. 

The period of the rain function can be set to match dlumal or sea- 
sonal intervals of high and low evapotranspiration potential. The 
circuit will act to increase the tln^ to soil saturation. or to raise 
the infiltration rate. Over this .long term, time variations in the 
infiltration rate can be ignored since the time required for it to 
fall to the final value is much shorter than the ev^otranspiration 
period. For the conplete model, a recurring rainfall with a period 
equ^ to the desired time between storms could be combined with the 
ev^KJtranspiration and infiltration circuits described above. . _ 

■ ^y- this -means, a rain event coiild' be 'initiate, ’then stopped arid a ' ‘ 
second event restarted. During the interim, the'ev^otrarispjtatlon 
• circuit would deplete soil nbisture;' 
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The analytic mcwlel is conposed of the follovflng modules: 

Rain recurrence module : This is an empirical fomalation, derived, 

fran analysis of the rainfall records of the 158 ARS test watersheds: 



Where: 

1 = rain rate, m/sec 
T = recurrence interval, years 
t = rain duration, hours 
ai, tt2 

03, i = constants, function of the location. 

Rain spatial correction module : For large watersheds, a spatial 

correction factor is introduced. This factor converts the point rain~ 
fall rate at the center of the watershed to a lesser effective rain- 
fall rate: 

K = cP 
e 

Where: 

Pq = effective rainfall rate, cm/hr 
P = point rainfall rate, cm/hr 
c = correction factor. 

The correction factor c is derived from the curves of Figure 30 . It 
shoiild be noted that the reliability of the spatial correction factor 
is ^ yet unproven for more than a few re^ons where data were gathered. 
It should thus be used with reservation. 



unvE 
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Rain Duration hrs. 


'• '■i * '• t 

' I 




Watershed 

Area 


Effective rainfall rate P© 
Pe= CP 

■ ..'i • 

P .* point rainfall rate 


40.000 ha 

20.000 ha 
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Subsurface abstraction module : Since initial infiltration rates are 

typically two to ten times the final rate occurring at the time of 
concentration and the processes are non-linear, accurate results 
require that the process be modeled over time, rather than relying 
upon some "average" value. The expression for subsurface abstrac- 
tion developed by Heggie Holtan of the AES has been arployed for the 
time being: 

• » 

I = GI • a • (Sa - + If 

\\ihere: 

0 

1 = total infiltration rate (or subsurface abstraction 

rate) 

61 = maturity of cover 

a = average vegetative cover factor 

Sa = average available water capacity = total available 
storage - initial moisture content 

I = cumulative infiltration 

If = final infiltration rate 

This formulation was selected after ccnparison with other existing 
formulations because: - 

1. Holtan’ s equation is being widely applied over a diverse 
spectrum of cover and soil types and conditions; thus, much 
in the way of oipirical results should be available for con- 
parison purposes. 

2 * It is a canplete formulation; it includes both surface 

(cover ccarposition and condition) and subsurface (soil type 



and antecedent moisture) phencmena. 


3 . Its results do not differ excessively from those of other 
widely used formulations, as shown in Section 3.2. 

4. It ejqjlicitly includes surface obserrvables (the a and GI fac- 
tors) which are potentially remotely sensible. 

• Overland and channel flov^ module ; As will be remaribered, the peak 
outflow from the watershed will occur at the time of concentration, 
assuming a rain of at least this duration and of constant intensity. 

At the time of concentration, all points of the watershed are contri- 
buting to the ruiK>ff . 

The initial analytic model for overland flow assumes the watershed 
schematic confi©u?ation shown in Figure 31. A single channel flows 
down the caiterllne of the watershed. Note that the single channel 
l^?pothesis is essoitially valid for subwatersheds of the highest 
order (the smallest in area); for larger watersheds, composed of 
several aggregates of subwatersheds, the contribution of each sub- 
watershed can be calculated with this model and the aggregate contri- 
bution is then con^iuted by routing. 

With reference to Figure 31 » consider the watershed subdivided into 
strips of unit width and extending perpendicular to, and fi?om the 
channel to the edge of the watershed. The total amount of water falling 
on such a unit strip must be released to the stream: 
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£* Average overlcmdflow 
length, meters 
N= friction coeff. 

S * Slope, meters/meter 
t « Rain durotlon, hrs, 

T» Recurrence Interval, 
years 

* R®Sional 

' xqlnfan 

co^ciems, 
derivobie 
from rainfoil 
records 

L* Channel length' jIjk 



IlMX' 


Qmox’ 


<=<> 

TtTdfs 

—JsfL, 

^06 §2/S 


I 


■ 2 tn r, 
% 


1 *^ 



^3URE 31 





Where: 


q = outflow, m^/sec 
i = rain rate, n/sec 
1 = flow length, m 


The flow velocity, which equals outflow divided ty cross section of 
flow, will be, TTiaintaining the assunptlon of unit Width: 

( 2 ) 


_ q _ il 

V = ^ = -j- 

a. d . 


Where: 


d = depth of flow, m 
V = velocity, m/sec i\', : 

a = cross section of flow, 


Equating this result to the Manning eqifetion, which also describes flow 
velocity, yields: ‘ 

. ( 3 ) 

■■ - -r- - ■ ■ 

Where: 

r = i^draulic radius • 

S = slope . (■ .-C‘. : - ■ 

n = average surface friction • factor 

I - . . * - 

i' i- 

Since, for overland (shallow) flow, the h ydra ulic radius equals the 
depth of flow, equation (3) becomes: 



' Whence; 


nil 
^ a (ISi) 

s'A 


Substituting expression (^) into the original equation for velocity- (3) 

• • . 1 ^ 

yields: 


•• 'J. 1:1 -^ i(il)^^5 


%J/io 


(5) 


Calling ’-'first time of concentration," T , the ’’ time required by the 
flo-w to traverse the entire unit strip of width 1; 

H 


T 


= i = (In) 
°(sec) V i^^s^° 


Converting time units to hours: 

3 / 




. = IiE) 


®(hrs) i^®s^/^“(3600) 


...X6X 


Substitutir:g the formula from the rainfall mbdOile for i in ec^ration 

r 

(6) yields: 


■'o- 


nUZ % 


(7) 


f^) s%"(3600) 


If (axlT'^^) is set equal to’^,"then: 


C 




1 


5 “(3600) 




(B) 


Replacing the t term in the formula from the rainfall module by 
equation (8), and assuming d to be small (as it is in most cases). 



•91“ 


yields the rain rate for the tinK of concentration duration: 


1 = 


T 



Or; 




I 

(36002 


1 


(9) 


5br a watershed vri.th a central charmel, such as that of Figure 31j 
the total overland flow length for the unit-width strip will equal 
2*1. In equation (7), then, 

q = 2il <10) 


Substituting equation (9) into (10) and summing along the channel 
length L gives: 


Q = 21U 


1 — (In)®*® “1 

5^‘*s^(3600T 


7 






( 11 ) 


Where: 

h = channel length, m 
1 = average overland flow length, m 


It is clear that this initial model assianes that the overland sur- 
face flow is much slower than the channel flow. This assunption 
turns out not to be overly in error, as the following more conplete 
version of the model shows. 


To add the effect of channel flow, consider the geometry depicted 
in figure 32. The total concentration time, which is now the sum 
of tte times required by surface flow and channel flow, is: 






Average overland flow length, m 

Nft= friction coefficient, ov&iond 
Ng® " “ , channel 

L ® length of chonner,m 
SM«,0verland, slopai m/m 

Sc»chonnd slo}^,m/m 
Tco» concentration tirhe, 
f^^AhannsI hydraulic radiusji'm 
d ® depth of overland flow , m 
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'total 




O CO 


a T 

^overland 


\1 + 


% 


K 


_0 

S_ Vz 
(~) 


(la) 


TOiere the sut>script “c" refers to the channel, and "o" refers to - 
the surface: 

R- * hydraulic radius, m 
c 

d * depth of overland surface flow, m 


The other parameters are as defined previously. The formulation (11) 
for the flow Q changes accordingly. 

The sensitivity of this model to observable phenomena is shown in 
Figure 33 for watersheds with various typical types of surface cover 
and varying drainage density. 


Point of flcxjding module ; In most practical applications, the user 
of a hydrologic planning model is interested not only in the accurate 
value of the peak flow, but in the coordinate viiere, along its length, 
the stream or channel will actually begin to flood. The stream 
begins to flood when the water level equals the hei^t of the banks. 
Flooding does not necessarily have to occur at the vxatershed’s 
outlet; it is a Jh^nction of channel shape, slope and rou^iness. 



0 0.01 


P|6U 

SURf 




151686 factors are conveniently combined into a single formilation 

e3^i*essed in t^mis of the most easily observable parameter, namely 

channel width, for a slnpllfied linearly increasing channel. 

Flooding begins to occur ^Aien channel width (w ) at any point along 

c 

the channel is less than: 

C O 

w. < U 

r+2k 

Where: 

L « channel length, m 
c 

? = roughness ratio surface/channel 

i 

k = channel geometric correction factor 
Yq = depth of overland flow, m 

w » channnel width at distance L frcan beginning of 
channel ° 

Figure 3 ^ supplies an -exanple for typical values of the parameters. 

The modules described above obviously cannot be singly connected. As 
is apparent almost from inspection, there are strong feedback factors 
between modules. The connection is thus best performed by program- 
ming the modules and their interconnections on a coir5>uting machine, 
analc^ or digital. 

Nevertheless, an analytical interconnection, using simplifying hypo- 
theses carefully checked for validity, was performed for a set of nine 
ARS test water^eds, chosen at random. The results of this ccstipari- 
son are presented -in the next section. 
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ON FLOOD GONDlTfONS 


Average Depth of Overland Flow- cm 






Lc=iaoocm 



3/8 e 3/8 ^ 5/8 


Wc< 


6 9 12 15 

Ctfflnnel Outfall Width-m 


\ I 4. 2k ' 

Lg" chanrel length-m 

S ' roughness ratio 
ov^iand/channel 

k - channel geonetric 
correction factor 

Yp- depth of overland 
flow-m 

Wg- channel width 
(5 outfall“m 

Example shown ; 
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5.0 DATA ANALYSIS AND MODEL VERIPICATION 

The modules of the .hydrologic .jjlanning model described in 
the previous section were interconnected analyticallya using' 
simplifying but carefully checked assumptions'. This was 
done on a randomly selected sample of the 158 ARS water- 
sheds, to gage the model’s applicability and accuracy and 
to 'test the hypothesis that valid predictions of peak flow 
events can be made -using remotely sensed data inputs. The. 
sample consisted of nine' watersheds distributed through- 
out the U.S. The output of the model, e.g. peak flow, 
was then compared: 1) to actual flow records and 2) 

to the predictions from the most popular current ungaged 
models . 

Three criteria were used to select the sample from the 
158 watersheds. First, a diversity of geographic and phy- 
sical characteristics was desired. . Watersheds were sel- 
ected 'from different areas of the country so that climate, 
vegetative cover, soil type, and other key factors gov- 
erning vfat'ershed behavior would vary as widely as pos- 
sible. Second, it was desired to include watersheds with 
the largest number of years of record. Thirdly and; fin- 
ally, watersheds were selected fpr which the most det- 
ailed’ topographic info rma-fc ion was available. 

A list of the 9 watersheds, ^their- locations and surface 
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areas follows: 




'x * *' ‘"J *' 

No. Assigned. 

Watershed . . 

Nomenclature 

Location 

Area km. 

1 

Watershed V/-1 

N, Danville, Vt. 

42.9 

2 

• - ' 1 * 
Watershed 19^ 

. ’ (* •* \ . * 
Coshocton, Ohio 

% 

3 

« • * »’’*r 

Thorne Creek 
’.Waiters he'd V7-1 

Blacksburg,- -Va.- 

* • ’. • . ji* 

M : ■ 

•Watershed W-10 

Oxfoyd, -Miss, 

22.3 

5 

Watershed V.^-l 

Pennimo r e , Wi sc . 

-1.3 

6 '■ 

VJatershed L21 - 

Chicas.ha, Okla., 

. 532.1 

•7 •. 

. Watershed D - 

Waco, Texas • - 

■ 4..5‘ 

8- ••• ^ ' 

Watershed W-1 . 

Safford,- Ariz. 

- 2.-1 

r 9 - it • - . 

Reynolds Creek • 
V/atershed W-1 

Reynolds, Idaho 

■ 233.2 


Average : 

Standard deviation a 

- 94.6 
: 180.1 


• obvious from the abov;e list that the selected, water- 

■ ■ • I • , • . b-.- '■ ‘ • V, ( „ 

sheds-,. vary xvidely in locatipn and area-. These. two fac- 

• tors alone haye rar’eat .influe.nce, upon watershed behavior: 
location determines vegetatiye. cover, soil typ.e, .bopo- 
-graphy.., and amount .of rainfall, area affects^ drainage den- 

•-.sity and -time of concentration,^. .. 


5.1 Analytical Procedure 




■ It should. e^ain -be emphasized that .the simplified analy- 
tical procedure hereinafter described is not intended to 


replace ulimately the"' c6’nstruc’t'ioh.'’o'f the ■‘compiete ‘model 
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■ coritlai'ni'rlg all modules properly ■'interconnected and integrated. 
THe ‘simplified analytical procedure was 'here used simply as 


a substitute for con^lete 'interconnect ion , which is sched- 
uled for subsequent phases of this project.. In spite of 
this ^shortcomings si^if“icant and most encouraging results 
were .achieved^, as will, .become, apparent from the. sect ions 
that .follow.... To determine capability to predict peak flow 
events from remote sensing data and to compare the new mod- 


el against existing models, the selected 9 watersheds were 
^aiyzed in detail by means 
'wiiich follow's . ■ ■" ^ 


‘■^aiyzed in detail by means of the step-by-step procedure 

^ ' id: i T O ; : i O • I . ^ i 


1. .. The. perc,en-tages pf,..each type of soil. -.within the 
.wa,t^r.sh,?d, . by ixifil^tratlpn class, were det.errained. 
..To th^s,_ end,,, Soil; Conservation. Service soil data 


were- used. SCS publishes a listing of the,. gener- 


al classification of all soils in the United States. 
'The SC'S cias'sifies soils into four general classes. 


designated A, B, 'C and D; 


/“i* J t . 




.Class , A... denotes .soils with- high infiltration rates 
evep...'VJh.e^..'fc}iQ??QUgbly;..w§|ite^ and, -therefore* with 
„.^ow runoff^, potential; ..Thpae. §oils ^ consist:. chiefly 
of deep, well-to-excessively drained sands .jor grav- 
els. 
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Class B denotes soils having moderate infiltration 
rates when thoroughly wetted'. They consist chiefly, 
of moderately deep to deep, and moderately well to 
well~drained, soils with moderately fine to. mod~ . 
erately coarse textures. 

Class C includes soils having low infiltration rates 
when thoroughly wetted. They consist chiefly of' 
soils vrith' a layer v;hich impedes downward movement. 
Of water, or soils with moderately fine to fine 
texture . 

Class D represents soils with very low infiltration 
rates vjhen thoroughly wetted, and, ‘ therefore', ' with 
high runoff potential. These soils consist chiefly 
of clay soils, soils with a permanent high water 
table, soils with a clay layer at or near' the sur- 
face layer and shallow soils over nearly imperv- 
ious material . 

An average watershed soil class was determined by 
computing a vreighted average of the above data. 

For example: for Coshocton, 86% of the soils are 

type C, while lH% are type B. Therefore, average 
soil class for Coshocton is approximately C, This 
average is used later in the conventional predic- 
tion f o rmul as . 
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2, An average infiltration rate was calculated for 

each watershed by means of the subsurface abstrac- 
tion module. This was done by first computing 

• • 

final infiltration rate (If) for each type of soil, 
then taking a v/eighted average . Values for final 
infiltration rates by soil class, presented in 
Table 11, were taken from the USDA HL 7^ model, 
authored by H. N. Holtan, G. J. Stiltner, W. H. 

Henson, N. C. Lopez, of the ARS (Reference .1 ) . 

The character of the soil vis-a-vis layers of soil 
which constitute an impediment to flow are used 
to determine the choice of the value within the 
range. A low value for impeding layer of clay; 
a mid value for loam; and a high value for sand. 

3* After average final Infiltration rate was calcul- 
ated, available water storage per unit depth was 
computed. Values of available storage capacity 
were assigned on the basis of soil type, according 
to^ Table 12. A weighted average was taken to deter- 
mine the available storage capacity (Sa). 

4. An average vegetative factor (a) was computed for 
the watershed. First, the distribution of cover 
■ was determined from the data base (for example, 

11^ cultivated, 58 ^ grassland, etc . . ). Then 

each type of cover was assigned a value according 

to Table 13- A weighted average was computed, according 







COVER GOCX?CONDITK)N PO(H CONDITION 

- : - 


Fallow = - ■ 

0 . 30 ^ 

0.10 

ROW Crops 

0^0 

p.io 

Small Grains " 

O.TO 

0.20 

Hay (legumes) 

040 

0.20 

Hi^ (sod) 

0.60 

‘ o.<o 

Posture ( bunch gross) 

0.40 

0.20 

Temporory . Pasture ( sod) 

0.60 

0 .40 

Permoient Posture (sod) 


- 0 -^ 

woods and Forests 

IJO 

0.80 
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to percentar.e of each type of vef;etatlve cover In 
the watershed. 

• 

5. The average infiltration I over time was then cal- 
culated from the equation: 

I = GI ’ a + If 

6. To enable comparison of results with the SCS procedure 
outlined in reference 2, an average SCS curve num- 
ber was computed. Curve number values were taken 
from Table 14 according to average soil group. 

(The determining factor in choosing a curve number 
is vegetative cover). A weighted average yielded 
a final curve number. 


■7 . The approximate time of concentration was computed 
applying the empirical equation developed by Kerby 
(reference 7) to each of the 9 watersheds: 

■2Ln1-'*^ 






Where : 


t(j = time of concentration 

L » distance from the most remote point in the 
basin to the channel, in a direction paral- 
lel to the slope 
S = slope 

n = retardance coefficient, according to Table 


15 . 
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Tabi@ 15 



Type of Surface 

■w ' 

Value of n 

Smooth im^vious surface 

0.02 

“ bare packet soil 

.0.10 

1 . Poor grass, cultivated row crops or | 

moderately rough bare surface 

.0.20 

Pasture or average grass 

,0.40 1 

Deciduous timbsriand 

0.60 

Conifer timberland, deciduous 
timberland with deep forest 

litter or dense grass . . 

V 

1 

0.80 



8, The surface friction coefficients (Manning's "n") 
were derived'Trom Table 9 (see Section 3) and a 
weig^|;ed average .calculated which included the ef- 
fects of the vegetati\)^e cover of the. ^watershed. 


9. A rainfall intensity-frequency-duratib'n relation 
was developed from NOAA published rainfall records 
■ by eiiipfrically fitting to each’ watershed the rain- 
fall module curve: ■ ' 





VJhere: 

i = rain rate, m/hr. 

T = recurrence interval, years . . 

t = rain duration, hrs . 

ai, ctz ' ' 

as, d = constants which vary vmth region. 

!{C: • • ' 'i- vi : -•.-.i •' 

As will be shown in Section 5.2;1, this formulation 
permits the determination- of* the"- intensity of a 
- 3 j -rain event .of -any -recurrence artd any -duration with 
. .:: 3 { -.;a high degree of acciaracy, to serve as the input 
into the hydrologic planning model,' 

10. Peak discharge rates per unit of watershed area 
. '.o • . • * . - - - 

were calculated from three .models in wide current 

usage, plus the new model,, in order to assess their 

•i*-&lat'ive accuracies.. - -These- ;mod.els, -were : . • . 
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a. Rational formula method (Reference 3) 

Q = CIA ' 

Where : ’ ' ’ 

Q = outflow, mvsec 

•• --c “ constant- based.-on soils and cover, dimen- 
sionless , % • ,,, . 

r i = rain, rate with .a .duration- o,f..T^,!and x- 

year recurrence, m/sec. 

A = watershed area, iii^ 

b. ‘ Cook’s method (Reference k) 

Q = f(A, R, I, C, S, P) 

V/here: 

f = an empirically derived function 
Q = outflow, m^/sec 
A = watershed area, ra^ 

•’■'••R =-wat’ershed relief ' factor, dimensionless 
I-= infiltration capacity factor,; ^dimension- 
dess .' = .q- 

C‘= vegetative -cover factor, dimensionless 
S- ^-'surface storage factor, dimensionless 
P = precipitation factor, dimensionless 

>;Sbll • Conservation Service' method' (Reference 2) 

Q '=^f(-A-;"c,''S,. T(j) ■ -V 

•-l^ere; - • . 


f = an empirically derived function 
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0 = outflov;, ni^/sec . 

2 

A = watershed area, m 

C = vegetative cover factor, dimensionless 
S = soil type factor, dimensionless 
time of concentration, sec 

d. Ecosystems method described previously. 

The results of these four methods were compared with the 
50-year recurrence flow (Q^q), derived from the statistics 
of the actual records using the Gumbel extreme value dis- 
tribution (Reference 5). 

! . 

5»2 Results of the Analysis 

The map. of Figure 35 shows the location of the 9 test water- 
sheds selected for detailed evaluation. Two watersheds 
among the nine - - Reynolds, Idaho (No. 9) and Chickasha, 
Oklahoma (No. 6) - - are very large, having areas greater 
than 200 km^. They are composed of numerous subwatersheds 
and,, therefore, do not precisely fit. within the framework 
of the simplified model described previously. As indicated 
in Section 3> watersheds composed of a significant number 
of subwatersheds require additional routing techniques to 
produce good accuracies. Their detailed analysis is res- 
erved for future phases of this, effort. 




figure 35 
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NINE WATERSHED S 
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5.2.1 Recurrence Rainfall 

For peak events, as previously described, the rain duration 
chosen was the watershed’s time of concentration. There~ 
fore, a rainfall intensity for a given recurrence period 
can be computed. 

The results of the empirical fits of NOAA rainfall data 
(Reference 6) to the rainfall module formulation for each 
of the nine v/atersheds are shown in Tables l6a through i-. 
Table 17 summarizes the results for the nine watersheds. 

The average error for the nine watersheds was 2.58^. This 
is well within the bounds of the 'errors of the other measure- 
ments. 

5.2.2 VJatershed Physiography, Vegetative Cover, and Soils 

The physiography for each watershed was developed from 
USGS 7.5 minute quadrangle, maps . It included the physical 
quantities: area, slope, channel length and drainage den- 

sity. The capability of remote sensing techniques to deter- 
mine these physical quantities will be covered in Section 
6, Table 18 summarizes the data for the nine watersheds 
shown in Figures 36 a-i. 

The vegetative cover and soil distributions were developed 
for each watershed as a weighted average using descriptive 
-data provided by the USDA-ARS. The utilization of remote 
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f?alnfall Daraflea-Rscurranc® Data- damviue, Vermont 


rOiiratlcnr.. 

Hre. 

- RECURRENCE , 

25yrs 50yrs 100 years 

* V > 

ae 

3in/hr 

.076m/hr 

3.4 in/hr 
.086m/hr 

3.5 in/hr . 
.089m/hr 

1.0 

'1 p 

I.SIn/hr 

,046m/hr 

2.lln/hr 

.053m/hr 

2.3in/hr 

.058m/hr 

r-t -s 

i? . 

2.0 . 

l.2ln/hr 

.030m/Nr 

l.3ln/hr 

.033rh/hr 

• ‘ ‘ i"> f ' ; 

1.4 in/hr 
.036m/hr 

• 3.0 

0.8ln/hr 

.020m/hr 

- ^J.;Qin/hr 

.025fn/hr 

f ,1 ln/hr_ J 
.028m/hr 

. 6.0 

0.5jn/hr 
.OI3m/hr =' 

'if'--. 

0.6in/hr 

.OI5rn/tr 

- ■! ; - ■ 

0.7in./hr 
' .OlSm/hr 

• j‘v ’ . 

lao 

0.31 in/hr 
..0079m/lf 

0.33in/hr 

.0084fn/hr 

0.^lnj4ir 

.0l07m/hr 

' '.‘T'V 

s^.o 

O.I9ln>iir 
,.004Sm/hr ^ 

• * ‘ *. 1 - ' 

0.21 In/hr 
_,-f - ,0p53m/hr . 

P; 7 _ I ■ 

' d.23jh/hr 
, . ..0058m/hr ,... . 


General Fomiula 
*■" (tt.2).77 


wh^e 

i= Intensity, /hr. 
T= recurrence, years 
t* durcjtion, hrs. 
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Tofel® Rainfall Dyratien- Rtetirren^ Data- coshocton, o. 


Durotion- 

Hrs. 

RECURRENCE 

25 yrs 50 yrs 100 years 

0.5 

3.4'Vhr. 

.086m/hr 

. 3.6 Vhr. 
.091 m/br- 

4.27hr. 

.I07m/hr. 

1-0 

• 

2j7hr. 

.053m/hr 

2r3'7hr. 

.058m/hr. 

2.67hr 
.066 m/hr 

2.0 

1.2‘Vhr. 
,031 m/hr 

l-47hr, 

.036m/hr. 

1.5'Vhr. 

.038m/hr. 

3.0 

f ' ' 

-93Vhr’ 
.024m /hr. ' 

I.07hr. 

.025m/hr. 

l.lVhr. 

028m/hr. 

6.0 

.557 hr 
.014 m/hr. 

.587hr. 
.015 m/hr. 

.637hr. 

.OI6m/hr. 

V \ 

.307hr. 
,0076m/hr I 

,337hr 
.0084m /hr. 

.387hr, 

.0097m/hr. 

1 

^ 1 

2^.0 

,177hr. 

.G043m/hr. 

1 

.I97hr. 

,0048m/hr. 

.2l7hr. 

.0053m/hr. 


®ral fom«to . . , 039 T- 15 where, ..... 

C t + 2) -83 i = intensify, m/hr. 


T= recurrence, yrs. 
t= -duration, hrs. 
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Table Rainfall Duration- Rseurrsnoa Data- ^‘-^cksbuto.w. 

I6e 


Durotion- 

Hrs. 

RECURRENCE 

25 yrs 50 yrs 1 00 years 

Q.6 

4.0 V hr. 
,I02m/hr. 

4.4 7hr. 
.1 12m/hr 

5.0V hr. 

.127 m/hr. . . 

1.0 

2.5Vhr. 
.064m /hr. 

2.9'Vhr. 

.074m/hr 

. 3.3"/hr, 
.084 m/hr 

2.0, 

I,57hr. ' 
,038 m/hr. 

I.SVhr. 

.046m/hr. 

2.0 7hr. 
.051 m/hr. 

3.0 

U3Vhr. 

,029m/hr. 

l.3"/hr. 

033m/hr. 

l.47hr. 
,036 m/hr. 

6.0 

.6?7hr. 

.0I7m/hr 

.75"/hr. 

,OI9m/hr 

.827 hr. 
.021 m/hr 

12.0 1 

.4l7hr. 

.OiOm/hr. 

.427hr. 

.OHm/hr 

.507hr, 

.0l3m/hr 

2^.0 

.21 VN-. 
,0053m/hr. 

.25Vhr. 

.C064m/hr. 

.277hr, 
.0069 m/hr 


,040 T''® 

{ t +. 2)-80 


i* intensity, m/hr. 
T» recurrence, yrs. 
t»duratjcn,hrs. 
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Tebl® Rainfall Duraflan-RseMrrsnos Data- oxforp. miss. 

I@d 


Durotion- 

Hrs. 

RECURRENCE 

23yrs . 50yrs lOOyeors 

0.5 

44in/hr. 
.112 m/hr 

5.0 in/hr 
.I27m/hr 

5.4ih/hr 

.|37m/hr 

1.0 

2.8 in/hr. 
.071 m/hr 

3.1 in/hr 
.079 m/hr 

3.4in/hr 

.086m/hr 

2.0 

l.75ln/hr 

,044m/hr 

1.85 in/hr 
.047m/hr 

2.15 in/hr. 
.055m/hr 

3.0 

l.27in/hr 

.032m/hr 

1.43 in/hr 
.036m/hr 

1,57 in/hr. 
.040 m/hr. 

1 

so 

. 

.77 in/hr 
.0l9m/hr 

.87in/hr 

.022m/hr, 

.97in/hr. 

.026m/hr 

leo 

^ in/hr. 

. .012m/hr 

.52 in/hr 
Oi3m/hr 

.57m/hr 

.OI4m/hr 

24.0 

.28 in/hr 
;0070m/hr 

.31 jn/hr 
.0078m/N* 

.33ln/hr 

.OOSSm/tv 


General Formula 

049 T’I^ where: 

j » (t»} ).75 ‘ ” intensity, m/hr 

T= recurr^ice, years 
t« diratlon, hrs. 
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Tabl® Rainfall Duration- Racarrsn^ Data- fennimore wise. 


Duration- 

Hrs. 

RECURRENCE 

25yrs 50yrs 100 years 

0.5 

% 

3.8 in/hr. 
.097m/hr. 

4.2 in/hr. 
.I07m/hr. 

4.6 in/hr. 
.I17m/hr. 

1-0 

2.4 In/hr, 
.061 m/hr 

2.71 n/hr. 
.069m/hr. 

2.8 in/hr. 

.071 m/hr 

: 

2.0 

l.4in/hr 

.036m/hr 

1,6 In/hr. 
.039m/hr 

1.8 In/hr. 
,044m/hr. 

3.0 

L03jn/hr 

.026m/hr 

l.i7in/hr. 
.0 30m/hr 

I.23in/hr 
.031 m/hr 

6.0 

.58 in/hr 
.015m/hr 

' 

.67in/hr ^ 
.OI7m/hr 

75 in/hr • 
.OI9m/hr 

lao 

.36 In/hr 
, 009m/hr 

.39 in/hr 
.OiOm/hr 

46 in/hr 
.OI2m/hf 

2^,0 

.21 in /hr 
.005m /hr. 

.23ln/hr 
•OOSm/hr . 

.26|n/hr 

,007m/hr 


. Gensral Formulo 

** + i= intensity, m/hr 

' * ' T- recurrence, years 

t= duration, hrs 
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Rainfall Dyratian-Raayrrtn^ Data- chickasha,okla; 


Duration- 


RECURRENCE 


Hrs. 





25yrs 

SOyrs 

100 years 

0.5 

6.0 in/hr 

5.8 in/hr 

6.4 in/hr 

.!27m/hr 

.I47m/hr 

.I63m/hr 

1.0 

3.2 in/hr 

3.6 In/hr 

4.1 in/hr 

.081 m/hr 

.09lm/hr 

,104 m/hr 

2.0 

1.85 In/hr 

2.2in/hr 

2.451 n/hr 

.047m/hr 

.056m/hr 

.062 m/hr 

3.0 

l.43in/hr 

l.57in/hr 

i 

IB in/hr 

.036 m/hr 

,040m/hr 

.G46m/hr 

6.0 

.83 in/hr 

.93ln/hr 

l.I in/hr 

.021 m/hr 

.024in/hr 

.027m/hr 

12.0 

.49in/hr 

.54in/hr 

.eSin/tr 

OI2m/tr 

.OI4m/hr 

.OI6m/hr 

2^^.0 

.283in/hr . 

.321 in/hr 

,358 in/hr 

.0072 m/hr 

.008! m/hr 

.0091 nn/hr 


General Formula i 

where: ; 

.17 i s intensity, m/hr 

i T= recurrence, years 

, -. 0 ^^ duration, hrs: 

(t +.2) 
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Rainfall Duration- Rasurrenoa Data- waco, texas 


Durotion- 

Hrs. 

RECURRENCE 

25yrs SOyrs 100 years 

0.5 

6.4 in/hr 
.I37tn/hr 

62 In/hr 
.157m/hr 

6.6 In/hr 
.I68m/hr 

1.0 

3.4in/hr 
.086 m/hr 

3.9 In/hr 
.099m/hr 

4.2in/hr 

.107m/hr 

2.0 

2.l0in/hr 

,053m/hr 

2.35in/hr 

.060m/hr 

2.65 in/hr 
.067m/hr 

3.0 

|.6in/hr 
.041 m/hr 

l.77iin/hr 

.045m/hr 

1.9 in/hr 
.048m/hr 

6.0 

.967irvhr 

.025nn/hr 

}.07in/hr 

.0Z7m/hr 

USin/hr. 

.030m/hr 

IZQ 

.567in/hr 

;0l4m/hr 

.642 in/hr 
.016 m/hr 

.73 in/hr 
.OI9m/hr 

a^.o 

.33 in/hr 
. ,008fn/hr 

.37lin/hr 

.009m/hr 

4l3m/hr 
.010 m/hr 


where- 

i *5 int^slty, m/hr. 
T« recurrence, years 
t = duration,hrs. 


Gena^l Formula 
-17 




.054 T 


IB 





SAFFORD, ARIZONA 


Duratlon- 

Hrs. 

RECURRENCE 

25yrs SOyrs 100 years 

0.5 

• 

2.8 in/hr 
.071 m/hr 

3.4in/hr 

097m/hr 

3.7in/br 

.094m/hr 

1.0 

L9 in/hr 
048 m/hr 

2.2 in/hr 
056 m/hr 

2.5 in/hr 
.064m/hr 

2.0 

I.!3in/hr 

.029m/hr 

* 

! . 2 5 in /hr 
•032m/hr 

l.40in/hr 

036m/hr 

30 

.80 in/ hr 
.020m/hr 

.92 in/hr 
.023m/hr 

l.03in/hr 

026m/hr 

6.0 

0.48m/hr 

OI2m/hr 

0.50in/hr 
.013 m/hr 

0.58 In/hr 
.01 5m/ hr 

12.0 

0.25 In/hr 
006 m/hr 

0.31 in/hr 
.OOSm/hr 

0.34 In/hr 
.009 m/hr 

2^0 

OI50in/hr 

.0038m/hr 

0.1 80 in/hr 
.0046m/hr 

.!88m/hr 

O048m/hr 


G^rai Rjrmuia: 

027T^® 
'• 5 ^ 


where 5 

I » intensity, m/hr. 
T= recurrerse, years 
t* duration. Its. 
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Rainfall Duratisn-Rssarrsnss Data- Reynolds, idaho 


Durotlcn- 

Hrs. 

RECURRENCE 

25yrs 50yrs .100 years 

0.5 

i.O In/hr 
.025m/1ir 

1.2 in/hr 
.030m/hr 

■■■■ " ■■ 

1.5 in/hr 
038 m/hr 

1.0 

.66in/hr 
017 m/hr 

.8 In/ hr 
O20in/hr 

O in/hr 
.023m/hr 

2.0 

.38in/hr 

.ODm/hr 

.47 In/hr 
.OI2m/hr 

.55 in/hr 
.OI4m/hr 

3.0 

.33in/hr 

.008m/hr 

.40 in/hr 
.OlOm/hr 

,43 in/hr 
.on m/hr 

6.0 

.21 In/hr 
.0053 m/hr 

.23in/hr 

.OO^m/hr 

.26 in/hr 
.0(^ nr/hr 

12.0 

.125tn/hr 
.0032 m/iw’ 

.I46in/hr 

.0037m/hr 

.167 In/hr 
.0042m/hr 


.075 in/hr 
.0019 fn/hr 

.O^in/hr 

.0022m/hr 

.100 in/hr 
.0025m/hr 


General Formula: 

.008 T*^ where: 

** iTri — i« intensity , m/ hr. 

T= recurrence, years 
t« duration, hrs 



r — 

ECO 

Na 

LOCATION 

<L I 

a 2 

<is 

d 

Average 
Error (%) 

m 

Danville, Vt, 

0.033 

0.16 

0.77 

0.2 

3,09 

2 

Coshocton, 0. 

0.039 

0.15 

0.83 

0.2 

3.34 

3 

Blacksburg, Va. 

0.040 

0.19 

0.80 

0.2 

1.8 0 

B 

Oxford, Miss. 

0.049 

ESI 

0.75 


1.37 

5 

Fenni more, Wise. 

0.04! 

0.14 

0.78 

— 

1.76 

6 

Chickasha,Okla 

0.055 


0.82 

02 

1.69 

B 

Waco, Tex. 

0.054 

0.17 

0.78 

OI 

3.31 

8 

Safford, Arl 2 . 

0.027 

0.19 

080 

0.1 

3.45 

9 

Reynolds, Ida. 

0.008 

0.25 

0.71 

0.1 

3.40 

Grand Average 2.58 


Oeneral Formula 


tti 


Where : ^ 

1 » rain rate, cm/hr 


T • recurrence period, years 
t« rain duration, hrs 
Gti, d2 


Average Error «ipj E e, where = 


ct3, d* empirical constants 

predicted “observed 


observed 





























































TABLE 18 PHYSIOGRAPHIC DATA SUMI^ARY 

FOR THE HIRE TEST WATERSHEDS 


ECO NO. 

LOCATION 

AREA km? 

— m:mi — 
AVG. SLOPE 

SHAPE 


m/nri2 

DRAINAGE 

DENSITY 

1 

Danville 

Vt. 

42.9 

0.12 

Triangle 

40,300 

1/1065 

2 

Coshocton 

Ohio 

076 

. .. 

0.172 

Square 

1,491 

1/510 

3 

Biackstxjrg 

Va. 

12.3 

0.123 

Ellipse 

L«2.6W 

26,900 

1/460 

4 

Oicfsrd 

Miss. 

22.3 

0.114 

Triangle 

20,273 

l/llOO 

5 

Fennknore 

WIs. 

■Bi 

0.00 

Ellipse 

L*2.IW 

2,364 

1/550 

6 

i 

Chiclocisha 

Gtla. 

532.1 

0058 

Ellipse 

L«2W 

76,014 

1/7000 

7 

Waco 

Texas 

4.5 

0.021 

Wedge 

L»2W 

6,618 

1/680 

8 

Sofford 

Ariz. 

2.1 

0.020 

Rectangle 
L« 5.3W 

6,195 

1/339 

9 

Reynolds 

Ohio 

233.2 

QI76 

Wedge 

L»1.9W 

164,177 

1/1420 
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FIGURE 

NORTH DANVILLE VERMONT ECO i 


Area 42.9 Km^ 

Slope J20 m/m 
Shape-Triangle 

Length of Channel 

40.3 Km^ 

Drainage Density 1/1065 
.16 - m/m^ 


APPROX true 

NorrH 


1 a 


.033 T* 


Ct + 2)'^^ 


Cover 

64% Cultivation 
I 7% Pasture 
\ 5% Idle land 
3% Homest^s . 

I % Rocxis 

Soils 

33% Stony silt loam 
1 8% Calais loam 
i 8% Royaiton loam 
I 6% Rocky loam 
7% Fine sandy loam 
4% Coirain fine sandy 
3% Peacham silt loam 
I % Misc. soils 
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R6URE 

36b 

COSHOCTON OHIO ECO-2 


Area ■ .76 

Slope • .172 m/m 

Shapes Square 

Lenqth of Channel - 1491 m 

Drainage Density* 1/510 m/m^ 

im/„= 03914 

(t <-.2)-83 


Cover 


Soils 


23% Hardwood Forest 
58% Grassland 
I i % Cultivated 
8% Miscellaneous 


33% Mu^lngum silt loam 
19% i^ene shallow loam 
17% Keene silt iKim 
17% Mixed silt loam 
14% Muskingum stony loam 
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THWE ora AT BLACKSBURG, VIRGINIA ECO- 2 

FIGURE 


2 

Area » 12.3 


Slope « .123 iJVfe 


^jspe « Ellipse; 1>2W 


Drainage Density » 1/^60 m/m^ 


Length of CSjanrKl = 26,6 tan 


i = 


.040T 


.19 


(t + .2) 


TSG 


Cover 6l^ Pasture 

32^ Cultivated (Com, Snail Grain, 
Hay Crops) 

H% Trees 
2% Idle Land 
155 Roads 



^ 2 


Soils 3055 Groceclose Silt Loam 
17% Lodi Loam 
14? Frederick Silt Loam 
12?.Lit2 Silt Loam 
■ 8? Greendale Silt Loam 
19 ? Other Silt Loams 


ECOSYSTEMS 
INTERNATIONAL INC. 
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d 


OXFORD MISS 
ECO 4 


Area 22.3 Km^ 

Slope .114 m/m 

Shape-Triangle 

Length of Channel 20.4 Km^ 

Drainage tensity l/HOO m/m^ 




.049 


.14 




Cover ‘ 

23% Cotton, corn 
B soy boons 
35% Pasture B 
idle land 
40% Woods 
2 % Bare gullies 


Soils 

50% Ruston Ind. 
sandy loam to clay loam 

16% Collins 

!4%Proy}dsm» silt loam 

l2%Loiing silt loom to 
silty cloy loam 

8% Grenada silt loam 





p 

Area »■ I.3 Kn 
Slope » .08 rt/a 


Shape = Pectanprie: ]>2.4i-J 
Draina^ Density = 1/550 m/ri^ 


Lenfth of Channel = 2.kH km 


i ca 



<t+.l) 


.78 


Cover 


23^ Com 
10^ Grain 
215S Kay 
23^ Pasture 
16SS Idle 
7% Roads 


Soils 50^ Tara Silt Loam 

19? Dubuaue Silt Loan 
23^ Dodrevllle Sdlt loar. 
8^ .Tudson Silt Loam 



A«»ROXrMAT£ 


TRUE NORTH 


lORIGIISrAL PAGE IS 
OP POOR QUALITY 


ECOSYSTEMS 
INTERNATIONAL INC. 






rjlICKflSiiA, m 6 




f 


N 

f 


Area = 532,1 Sti 
Slope *= .058 nv'm 
Shape = Ellipse; L=2W 
Lenfrth of CJiannel = 7^,4 kra 

DrainaA^ Density = l/JOOO nV'm' 
i ^ ...055T;^ 


tftmxwtytjm 

DBCUtonoN, 


(t + .2) 


.82 



Cover 8% Alfalfa 

48? Sowed Crops 
44? Row Crops 


So ils 40? Noble Cobb Loam 

36? Darnel Vfoociward Quinlan Fine y 

Sandy Loam 

16? Port-Yahola Pulasici Fine Sandy Loam 
8% Noble Vanoss Cobb Sandy Loam 




FcosrsTr 

lNTr-:RNATIONA( 
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FI6UR 




RIESEL(WACO) TEXAS ECO- 7 


Area 4.5 Km^ 

Slope .021 m/m . 

Shape W L= 2W 

Length of Channel 6.6 Km^ 

Drainage Density 1/680 mAn^ 

(t+i) 78 
Cover 

60% Pasture 
6% Small grain 
. 3% Corn 
7% Cotton 
9% Row grains 
2% Gravel 8 paved roads 
1 3% Otter, mostly weeds 

Soils 

66% Wilson clay loam 
24% Burleson Heiden clay 
4% Frio clay loam 
3% Crockett loam 
2% Burles(^ clay 
i % Houston Black clay 
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PteURE 36h 

SAFFORD, ARIZONA ECO-8 


Area 2.1 Km2 

Slope .020 m/m 
Shape Rectangle L« 5.3 W 
Length of Channel 6.2 K 
Drainage Density 1/339 m/m^ 

. ( tV .|).8 


Cover 

85% Bore 
1 5 % Grosses 


Soils 

47% Tragus 
4% Gilman 
3% Laveen 
46% Luzena 
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FI6URE 36i 

REYNOLDS CREEK IDAHO ECO -9 


Arw 233.1 Km^ 

SiSpe .176 m/m 
Shape- Wedge La|.9W 
Length of cht^el^ 164,6 Km . 

Drainage Density - 
;^.008T-^® 

Cover 

95% Sagebrush & Rangeland 
2% Small stands forest ' 
3% Alfolfa - . . 


Soils 

12.65% Reywat 

1 0.68% Narmeli 
8.73% Bak«)ven 
633% Gabeca 
7.62% Ruclick 

6.87% Takeuchi 
5.04% Nannyton 
3.75% Lasing 
3.7i%Gemid 
3.08% Babbingtoh ’- 
3.01% Searia 
2.95% Glasgow 
2.90% Farrot"'''^' '^' ’'''^ 
2.81% Kanlee 
2.19% Castle Valley 
2.14% Nettletoo 

I 3.55% Additional 
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sensing for the development of this information is discussed 
in Section 6. The data for each of the basins analyzed is 
presented in Table 19. 

Using the above data, the pertinent parameters for the four 
models to be compared were computed by the procedure pres- 
ented in section 5.1. The parameters computed were: 

1) the time of concentration; 

2) the average infiltration; 

3) the average surface friction; and 
the average subsurface abstraction 

( 

Additionally, some of the four models, specifically the 
SCS model, required the computation of other parameters. 

A summary of the results of these computations is given 
in Table 20. 

5.3 Results of the Model and Comparison with Actual Rec- 
ords and Existing Planning Models 

Having thus determined analytically the values charact^er-' 
istic of the recurring rainfall, and having approximated 
the subsurface abstraction, the difference between the two 

t 

became available as the excess rainfall,, which is the gen- 
erator of the overland flow. The overland flow itself 
and the corresponding peak flow was then computed by means’ 
of the overland flow module detailed in Section 4. 
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table 19 


, ECO NOV 
i LOCATION 

Soil Class % 

1 A 1 

B 

C 

D 

I. Danville Vt. 

B 

45 

18 

36 

2, Coshocton Ohio 

0 

14 

86 

0 

3. Blacksburg Vg. 

B 

39 

61 

0. 

4. Oxfo'd Miss. 

0 

50 

50 

0 

5. Fennimore Wise. 

0 

100 

0 

0 

dChickosha, OKia. 

B 

100 

0 

m 

7 Waco Texas 

0 

B 

0 


8. Safford Ariz. 

B 

B 

0 

93 



; 1 

1 




9. Reynolds Ohio 


0 


17 


60 23 


SURFACE COVER 


@4% Hsrdscod fonisf}i6%pa8ft!79>«{Bli 
17 % Cultlvotod fK»yt 1% row crept 
9 % Dtmw bru^ oi^ dress 

23% Ka'dnood forsefy 88%6ra83toitd} 

I! % CdltlvQted} 3% mite. 


89% p^tursi 3% corn, 26%0(iltivored {toy, 
4% forestt S% (dl9 (^ea^t 1% imporis^e^ 


29% Rowereps, 39% pasture 3 fotlowt 40% 
forest, 2% base 


23% Corn, 10% drain, 21%(ny, 23%posture, 
I6%fal!cw, 7%ro«ls 


^% erelRcropt, 

44% Rowoit^ 

3%Alf0ffo . 

60% Pasture, 6% dre^n, 3%oore, 7% ootfoiH 
9%Retsorop9, 2%roedt, i9%Wseds 


36% Sort, 16% gross 


99% future, £%fer«9t, 3%Ai?e{ffa 





















































TABLE 20 SUMMARY OF COMPUTED FLANNIK© l«ODSL DATA 

FOR THE TEST Y^ATERSHEDS 


ECO NO. 

LOCATION 

TIME OF 

CONC. hrs 
( i ) 

SUBSURFACE ABSTRACTIONS | 

CN (4) 

0/P 

I^Soce 

Friction Coef 

Iq cm/hr. 

« 

If cm/hr 

(2) Sa m/m 

(3a (Holton) 

1 

Danville, Vt. 

1.6 

8.3 

0.33 

.24 

082 

. 

71 

046 

0.079 

2 

go.^octon. 

039 

5.3 

0.30 

.20 

0.56 

74 

037 

0.048 

3 

Blacksburg, 

Va. 

0.50 

7.2 

0.38 

.21 

072 

75 

Oil 

0.039 

4 

Oxford, 

Miss. 

17 

6.4 

0.41 

.23 

0.56 

71 

0.18 

0060 

5 

Fenninnore, 

Wise. 

0.53 

5.0 

0.58 

.23 

0.42 

70 

0.14 

0034 

6 

Chic kasha, 
OKIa. 

1.9 

2.44 

048 

.26 

016 

76 

0.04 

0.038 

7 

Waco, Tex. 

0.72 

6.10 

0.08 

.21 

0.65 

84 

0.18 

0.035 

8 

Sofford, 

Ariz. 

0.9 

1 

2.06 

0.08 1 

.20 

0.23 

94 

005 

0.020 

9 

Reynolds, 

Ohio 

1.28 

1 

2,95 

0.23 

.20 

0.31 

74 

0.|9 

0.036 


(i) BY KERBV'S EOUATIOM (2) AVERAGE Y/ATER STORAGE CAPACITY (3) EMPIRICAL C0W8TANT I M HOLTANS IKFILTRATION 

FORMULA 


(4) 80S CURVE m, CK 
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The peak flow thus predicted from the new model - - which 
for ease of recall we shall hereinafter refer to as the 
ECO model - - was then compared .to each of the nine water- 
shed’s streamgage records. 

Finally, the predictions of the other three planning mod- 
els in widest current use - - the Rational formula method. 
Cook’s method, and the SCS method - - were computed for 
the nine watersheds under the same conditions, and their 
results compared with the streamgage records and with the 
predictions from the ECO model. 

5.3.1 50-Y'-r.r Peak Plow from ’Actual Records 


The mean peak flow and its variance were developed from 
the records for each test watershed. The Q^q, i.e. the 
peak flow corresponding to a 50-year return period, was 
calculated using Gumbel’s extreme frequency distribution: 




Where: 

Q50 = peak flow with 50-year return period 
^ = average yearly peak flow 

e variance of the average peak flov; 

K = Gumbel constant which is a function of recur- 
rence period and length of - record. 


The values of K are presented in Table 21. 


TABLE 21 

VALUES OF K FOR DIFFERENT SAKAPLE SIZE B RECURRENCE 


Sompie 

size 



Recurrence Intervol 


. , 

10 

20 • 

V 25 

50 

75 

100 

is: 

1.703 

2.410 


3.321 

3.721 

4.005 

20 

1.626 

2.302 


3.179 

3.563 

3836 

25 

1.575 

2.235 

2.444 

3.008 

3.463 

3.729 

30 

1.541 

2.188 

2.393 

3;026 

3.393 

3653 

40 

1.495:; 

2.126 

2.326 

2.943 

3.301 . 

3.564 

,.50 

1.466 

2.006 .. 

2.283 . 

2,889 

3.241 

3491 


1.446 ^ 

2.069 

2.253 

2.852 

3.200 

3.446 

; -70 

1.430 

2.038 

2.230 

2.824 ' 

3.169 

3413 

75 

V 

1.423, 

2.029 . 

2.220 

1 » 

2:812 ^ 

3.155 

3.400 

100 

1.401 ' 

1.998 

1 

2.187 

r , 

2.770 

1 

3.109 

3.349 
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For example, to compute the 50-year event from a data base 
of 30 years of peak flow data, the corresponding value of 
K from the table would be 3.026. Therefore, in this case; 

*^50 “ Qp 3.026 UQp 

Table 22 summarizes the results of the statistical analysis 
of the records. It was used as a basis to compare the ECO 
model and the three planning models in current practice. 

5«3.2 Remote Sensing Model (ECO Model) and Results 

The ground rules and assumptions employed in applying the 

f 

simplified analytic formulation were: 

1. The expression for the time of ooncentrat ion must 
be made to represent the maximum time required for 
precipitation to reach the outlet from its point 
of impact. 

2. Infiltration abstraction can be accounted for, by 
modifying the rainfall accordingly . 

Referring back to the procedure detailed in Section 4, the 
time of concentration can be described as the maximum of; 



V V 

o c 


( 1 ) 



TABLE 22 Q 5 Q 


Fl«)M STATISTICAL ANALYSIS OF RECORDS 



ECO 

NO. 

LOCATION 

AVERAQE YRCr. 
PEAK FLOW 
5p, 

PEAK FLOW 

VflRIAKCE 

0"Qp,«?>&ats/Iair 

YEARS OF 
RECORD 

Oao 

®90 

mVMoAn^ 

1 

DANVILLE ,VTi 

.38 

,169 

10 

3.04 

.98 

2 

COSHOCTON, 0 

2.46 

2.28 . 

8 

3.60 

10.8 

9 

QLACK8BURS,VA. 

.32 

■ • .30' 

III 

3.40 

1.3 

4 

OXFORD, MISS’ 

3.18 

281 

II 

3.^ 

12.0 

S'*' 

FENNIMORE, WI8C. 

2.88 

/ ♦ 

2.98 

20 

9.09 

11,8 

6 

CHiCKASHA, OKLA. ' 

.17 

» 

.20 

9 

;^80 

89 

r 

HmCOiTEX. 

3.80 

9.29 

£0 

3.09 

13.8 

0 

jBAFFORD,ARIZ. 

1.89 

1.48 

29 

3.04 

9.2 

—J 

9 

REVNOtDS, IDAHO 

.18 

49 

9 

1 

1 

989 

■ .87 


I 

H 

<jO 

00 

1 
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I'Jhere: 


= overland surface velocity and length 
V , 1 = channel flow velocity dhd length 

w C 


Prom Manning’s equation: 


VJhere : 


. _ ^o ^o , '^c ^c 

o o c c 


(2) 


n , d s = overland surface friction, depth 
o o o Qj. flow, and slope 


• ^c» ®c ” channel friction^ .depth -o'f 

and slope 


-frowp' 


Expressing channel characteristics as a function of sur- 
face 'characteristics: 


■■ "0= Vo...... 

*c‘ =' Vo “ ' 

. 

=0 = ‘"4=b 




( 3 ) 


V/here k^, kg, k^, k|j are constants, 


Substitution of these terms into equation (2.) yields: 


io". 


^ , VoV o 


-total 


• • ■ • i-f'-yxo ‘ : 

O 0 




. (^0 


f'r. L! 


’OESeiHMJ PAGE IS 
OP POOE ^JirALFTYi 



= t. 


'total 


1 + 




Or: 


'total 


Where: 


(J> t. 




( 5 ) 


/ + 2 / 1 / . 

The values of k^j k 2 , computed from flov; 

records and from knowledge of watershed characteristics. 

The correction for subsurface abstraction -is made, by re- 
ducing the rainfall input by the amount of water which 
becomes infiltrated. This was accomplished by introducing 
a factor into the rainfall expression to degrade the 
rainfall rate, leading to the modified expression: 


K T' 
i = -2^ 


a 2 






as 


Where; 

k - % of, -rainfall which does not infiltrate 
i = pi^ec^plt'atdon excess 


It is worth emphasizing again. that this expression is 
only approximate, sinc'e in reality- the relationship be- 


tween rainfall and subsurface abstraction is non-linear. 
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Thus, the results expected are necessarily approximate; 
more accurate results should be available from the com- 
plete model. 

The k term above was calculated by -comparing the rain 
rate occurring over the time of co.ncentration with the 
average infiltration rate for the same period. For 
example, the rate (P) for the 50-year recurrence, t 

<2 

duration rain within the Blacksburg watershed is 0,109 
meters/hour. The infiltration equation, using the con- 
stants for this watershed, derived as explained in Section 
5.1, is: 

^cm/hr " *72(4.97 - + .38 (7) 

. ' ■ -• a. - 

For t = , 52 hours (from the Kirplch formula), the infil- 
0 

tration rate will fall from 7.1 cm/hr to 3.4 cm/hr, 
with, an average value equal to approximately ^.2 cm/hr. 

The k factor, therefore, for this case equals: 

‘ “ i‘ 0.0i»2 m/hr 

1 - . = 1 - = 0.62 . ( 8 ) 

P 0.109 m/hr 

In other words, for’ this particular rain event, and for 
the Blacksburg watershed, approximately 62? of the rain- 
fall becomes runoff. 

When both ground rules and assumptions are included, the 
ECO formula for the overland flow module becomes: 
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• Q « 2L1C 


<P (In) 

s"/i 0(3600) 


1 

Vs n^.n-%3 


(9) 


Vlhere: 
% = 


k T 
ai 


Ct2 


Table 23 supplies values for the constants developed above 
for the nine test watersheds. The 50-year recurrence 
flows, as calculated from the ECO formula, are presented 
In the final column. 

5.3.3 Comparison of the Four Planning Models 

The predictions of the ECO remote sensing model in its 
simplified form were compared to the predictions of the 
three principal planning models in present use for un- 
gaged watersheds. The predictions were calculated using 
exactly the same base parameters for all models, namely 
the ARS and USGS data as described in Section 5.1. Com- 
puted values for the overall results are tabulated In 
Table 24, and depicted in Figure 37 for the nine test 
watersheds , 

In Figure 37, the solid line of unit slope (45®) is 
the locus of the points for which the predictions equal 
the measurements. Points falling above this line are 

f 

underestimates; i.e., the predictions fall short. Points 




o o 

sg 

wg 

gg 

’g 


ECO 

NO. 

LOCATION 

m 

MAlfl CHAN- 
NEL LENGTH 

LtRl 

OVERLAND 
FLOW LENGTH 
l,ni. 

0 

n 

i 

®80 

nS/aec/itwd 

1 

DANVILLE, VT. 

n 

93H 

eos 

2.0 

.079 

.12 

.91 

& 

C03K0C1X){>I,0. 

l.4*ICf® 

SIS 

■ 

163 

1.73 

.048 

.172 

26.S 

3 

SLACK39UE^,VA. 

■ 

SOOO 

167 

I2.I 

.03® 

.123 

1.01 

4 

OXFGND,£3IG3. 

I.4kKT® 

12,000 

€S0 

2.61 

.060 

.114 

10.8 

e 

FENNI^ClStE, WI3C. 

1.43 10**® 

1045 


3,23 

.0254 

.080 

12.8 

G 

ailCKASHA,OKLA. 

a.osjQi^ 



34719 

@00 

34.4 

.033 

.008 

.08 

7 

WACO TEK. 

1 

! 

■ 

3290 

263 

3.1 

.030 

.021 

H.5 

8 . 

1 

! 

8AFFODD,ARIZ. 

|.2kI0“* 


133 

6.3 

.020 

i 

,020 

8.3 

9 

! 

NEVN0LD8 t0A»0 

■ 

2HSI 

189 

32.2 

,038 

,176 

.001 









































\, Danville Vt. 

2. Coshocton, Oh© 

3. Blacksburg -Va 

4. Oxford Miss. 

5. Fennimore Wia;, 

6. Chickasha, Okia 
7 Waco, Texos 

8. Safford, Ariz. 

9. Reynolds, Ohio 






^ CKICKA8HA • eowpottd of tabbotlM 
Oato o« 9uMa*)n» net evoilofrio 
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falling below the unit slope line are overestimates; 
the predictions are too high. The dashed lines around 
the unit slope line indicate the region bounded by + 15/S 
error around the measured, or “true," quantities. 

It can be seen that in most cases the ECO model, even- 
though only used in its simplified form, is an improve- 
ment over conventional models. 66^ (6 out of 9) of 
its predictions are contained -within the + 15/S error 
bound. Two of the three predictions which fall outside 
the + 15^ bound, specifically Chickasha, Oklahoma, 
and Reynolds, Idaho, pertain in reality to ve 2 ?y large, 
complex watersheds which require a somewhat different 
procedure. Namely, that the model be applied to each 
of the individual subwatersheds of which they are com- 
prised, and that the outputs be then coalesced by routing. 
This approach is reserved for future effort. 

It would of course be premature to claim that the results 
from this limited sample provide proof that the ECO model 
concept is valid for all regions of the U.S. and for all 
flood regimes. Rather, the results tend- to support the 
validity of, and encourage the approach of constructing a 
hydrologic planning model highly sensitive to remote sen- 
sing data inputs. Further effort is required to deter- 
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mine the sensitivity of this type of model to its key 
parameters: drainage density, variation of slope, basin 

area, multiple watersheds, etc. The approach does, how- 
ever, appear to offer promise of yielding a practical 
model capable of using satellite data inputs, particular- 
ly in the future complete version. 

The following section describes the relationship of 

conventional and future reirate sensing techniques to 

/ 

the hydrologic planning model described. 



INTENTIONALLY 

LEFT 


BLANK 
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6.0 THE RELATIONSHIP OF REjyPTE SENSING TECHNIQUES TO 
HYDROLOGIC FLAMMING IVPDKLS 

The preceding analysis indicates that the surface features of a 
watershed play a significant role in the prediction of the peak 
runoff. This section addresses the applicability and feasibility 
of measuring surface features, and infeirring subsurface character- 
istics, by remote sensing techniques. 

Virtually all the requirements for the measurement of the geometry 
of surface features have already been adequately matched by remote 
sensing from aircraft. It remains to assess whether spacebome 

j 

sensors of modest geometric resolution but with high radiometric 
content limit the accuracy of the data. 

The requirements for the identification of the types of species 
on the surface hinge upon multispectral techniques of discrim- 
ination. Techniques to satisfy them are currently in the advanced 
development stages. What needs to be assessed is whether the 
identification accxrracies currently being esqperienced are consonant 
to the accuracy of prediction required from hydrologic planning 
models. 

Table 25 depicts the Inportant information required for the develop- 
ment of hydrologic planning models onploylng remote serving 
techniques . 
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6«1 Rgnote Sensing Observations 

Several of the information requiranents of Table 25, such as 
ai*ea and channel width, can be met by direct ji^asurement from 
ranote sensing Imagery of digital tapes. Other information 
requiranents result from combining two or more remote sensing 
measurements. For example, drainage density is measured by 
using directly observed drainage pattern plus area mensuration. 
Table 26 summarizes, for each element of information, the existing 
technique(s), the degree of feaslblity already demonstrated, the 
principal researchers in each technique, and the typical numerical 
values of accuracy attained at present frcm processing of 
ERIS data. 

6.2 Relationship of Remote Sensing Observables to Hydrologic 
Planning .Model Requirements 

The modules developed In the previous sections fundamentally 
require six elements of data: • watershed area, overland surface 
friction coefficient, drainage density, channel width, subsurface 
abstraction and slope. Significant information is currently 
available, in the U.S. and developed nations, on four of these 
elements: area, surface friction coefficient, drainage density, 

channel width. Additionally, a considerable boc^ of direct and 
inferential evidence exists frcan which to deduce the subsurface 
abstraction ccnponent. The last element, slope, can be addressed 
in -three ways: 
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1. ^'Jhere topographic maps exist, overlays can be made and 
adjusted to known benchmarks. 

2. Use can be nade of the existing Defense IVfepping Agency 
topographic digital tapes (which cover the entire U,S. 
at 100 meter contours). 

’ 3* Where there is overlap in the ERTS pictures, stereo 

pairs can be developed from which to measure the slope. 

In the future, this stereo capability may become 
routinely available in advanced Earth Observation 
Satellites - 

Table 27 synopsizes the techniques to extract the information 
required by hydrologic planning models fixim the remotely sensed 
observables corre^onding to each information element. 

6.3 Visual Interpretation of an ERIS Image of a Test Watershed 

A test was made to determine the extent of the information 
requirements of Table 25 which can be satisfi^' by analysis of 'Ms 
Imagery, figure 38 shows a ^x enlargement of a section of a 9** 
Band-S" ERTS transparency. Shown is the test watershed at CMckasha, 
Oklahoma, and its surroundings, figure 39 gives a corparable 
USGS topographic msp. The magnification was chosen to match' 
the scales of the two maps (1:250,000). 
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RELATIONSHIP BETWEEN REMOTE SENSINa OBSERVABLES 
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Observable 
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Technique 
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Drainage 
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coefficient. 
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Channel width; 

Measure by using "pixel splitting" in high' contrast 
situation between water and surrounding land surface. 

Slope 

Apparent relief. " 

Stereophotpgramraetry techniques on overlapped imagery. 

Channel Plow 
Capacity 

Channel width. 
Drainage density. 
Drainage pattern. 

Use existing enplrical relationships between channel 
width, channel meander, drainage pattern to esti- ; 
mate bankfull capacity of channel as a function of 
channel length. ■’ 

Soil Pennea- 
bility and 
MDisture ^ 
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Land use. 

Use soil associations - by classification. 
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TOie general outline of the watershed is apparent in the upper 
left and central areas of the ERTS image. The delineation of the 
exact boundaries, of the basin was acconplished in this case by 
conparison with "the topographic inapj simply because the map was 
.available. However, where maps are unavailable or of doubtful 
reliability, two. methods’ exist, for determining watershed boun- 
daries directly iVom the ERTS image: 

1. The segnents lying between streams which drain into the 
watershed undo? study and those which drain away from it 
into other water^eds may be divided in some ^yeighted- 
fashion. If no other infoniiation Is available, division 
can occur at their mid-point. Referring to the ERGS 
Image, streams which drain into neighboring basins can 
be seen at .points A throu^ I, while the streams of the 
test watershed are shown in tte overlay. By this means, 
sev^al reference points are established and can be 
connected to fom an approximation of the watershed area. 

2. In regions where the land relief is pronounced, ridge lines 
are visible and can be followed to delimit the drainage 
area. The average slope in the Chiclasha watershed is 
less than 6^, however, so this method is not reliable for 


this case. 
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The detennination of drainage pattern and density builds upon 
the information derived above. Referring to the Image and the 
overlay, the central channel is visible as a light gray strip 
running almost North to South. Secondary and tertiary channels are 
also visible. 

An advantage of EE?TS imagery is apparent here. This image was taken 
in October of the year when vegetation density is low. This 
makes obvious streams which might not have been apparent in the 
aerial yiiotographs from which the topographic maps are made. 

This appears to be the case at this location. More streams can 
be seen fran the ERES image than are recorded on the map, 
yielding an improved measure of drainage density. Stream length 
may be measured directly from the EKES image and divided by 
watershed area to yield drainage density. 

Drainage pattern, in this case one central channel branching into 
several ancillary streams, is immediately apparent. A second 
advantage of ERTS Image analysis is clear. The USGS has noted that 
a majority of the secondary streams of this watershed are ephemeral, 
that is, water does not always flow in than. Actual drainage 
density, therefore, may change seasonally. Consequently, unneces- 
sary errors could be introduced into hydrologic models. The high 
frequency retuam period of EESS provides the cs^jability of measiire- 
ment of drainage density at time intervals adequate to insure accuracy. 
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An litportant measurement is the determination of channel width, 
which,.* as noted earlier, is a detenninant of flood potential. 

The channel boundaries are visible on the 4x enlargonent and 
more apparent on the 8x magnification (Figure 40), A precise 
measuranent would require sub-pixel processing or a quantitative 
examination of the contrast between land and water on the image. 

A combination of drainage density and pattern and of channel width 
information will srleld an Indirect measure of channel carrying 
capacity, as described in Section 4. 

To measure surface friction ERTS imagery can improve surface cover 

i 

estimation. For exanple, ABS records -for this watershed show 
cover to be 48^ sowed crops, 44^ rovr crops, and B% alfalfa: 

This most recently published data is frcra 1967. It is clear fran 
the image that changes have taken place since; for example, 
about 20% of the watershed is covered with forest, exaiplified by 
the dark gray areas surrounding streams. This is confirmed by the 
location of shaded (forest) areas on the USGS map. Further, changes 
in cover can be expected to take place seasonally. These fluct- 
uations can cause up to a tenfold change in overland flow velocity. 

Other types of cover are equally visible. Surface water appears 
as very dark areas on the positive. Large bodies such as Fort 
Cobb Reservoir (lower left) and Lake Chlckasha (lower right) are 
exanples. Smaller bodies form black spots, as called out on the 
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overlay. Here again, it is clear that sane changes have taken place 
since the topographic map was assembled. Pbr exanple, a new 
lnpouncMent exists on Wildcat Creek (point J). . 

Urban areas can also be noted. Binger,, a town of population 603, 
lies within the watershed and can be distinguished as a rectang- 
ular area at point K. Even more geonetric regularity is present 
in cultivated areas. 

No effort was made at this juncture to infer subsurface para- 
meters from the EFEPS image. The subsurface parameters are in 
this case available from records. In combination vjith such infor- 
mation, it appears that techniques for extracting information from 
ERTS imagery, when fully operational, could play a significant role 
in supplying data required by hydjxi logic planning models. 
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7.0 CONCLUSIONS 

1. - An improved model for the prediction of peak flow 

events has been structured, which is specifically 
designed to take maximum advantage of the data and 
information stream available from remote sensing. 

2. The development of the model has been carried to the 
point where the overall framework has been con- 
structed and .five modules simulating the behavior • 

of significant hydrologic processes have been developed. 

3. The improved model is considerably more sophis- 
ticated than conventional hydrologic planning models. 

In particular, its modules are not simply inter- 
connected, but require feedback. In spite of this 
greater complexity, however, the model- is readily 
adaptable to analog computation with modest amounts 
of hardware. Preliminary sizing shows that the 
technique can also be programmed onto one of the 
smaller types of digital minicomputers. 

4. The model was exercised — not in its fully inter- 
connected form, but rather in a simplified version — 
to predict the peak runoff from nine experimental 
Agricultural Research Service watersheds, selected 
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at random from among a set of 158 instrumented and 
well -described watersheds. 

5 . The predictions of the new model in its- simplified 
version were tested against: 

a. The predictions from three of the most employed 
contemporary planning models — i.e., the Rational 
formula method. Cook's method, and the Soil 
Conservation Service method, 

b. The streamgage records of the nine test water- 
sheds , 

6. The results indicate that, within the range of 
applicability of its simplified version, the new 
model appears to be considerably more accurate 
than conventional hydrologic planning models. 
Specifically, in six out of nine of the watersheds 
tested, the new model supplied predictions of peak 
flow for the 50-year event falling within error 
bounds' of + 15^. For these same six watersheds, 
conventional models yielded discrepancies with res- 
pect to the records ranging from a minimum of 

1,2 to 1 to a maximum of 15 to 1. For the 3 remain- 
ing watersheds, the new model yielded predictions 
of lesser accuracy — the worst being 2 to 1, 
Reasonable explanations for the discrepancy are: 



a. The fact- of having oversimplified the model by 
not operatiri£5 it in its fully Interconnected 

-o 

version. 

b. The three watersheds are considerably more com- 
plex than the other six, and they need' to be 
split into subwatersheds, predicting? the output 
from these, then routing; all outputs through 
the watershed channels. This technique, which 
appears to be -well in hand, is proposed for 
future phases of the effort. 

The appropriate techniques whereby to extract the 
inputs and parameters required by the new model 
from remotely sensed information — whether imagery 
or digital tapes — were explicitly defined. Their 
feasibility was identified from specific past and 
ongoing ERTS investigator efforts. 



8.0 APPENDIX 


The Appendix, which summarizes pertinent data for the 
158 ARS test watersheds, includes the following infor- 
mation: 

® Watershed number 
° Location 
® Area 
® Slope 
® Shape 

® Shape correction factor 

® Time of cone ei>trat ion in hours and minutes 
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